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We construct a map between the quantum field theory of free Weyl or Majorana fermions and 
the probability distribution of a classical statistical ensemble for Ising spins or discrete bits. More 
precisely, a Grassmann functional integral based on a real Grassmann algebra specifies the time 
evolution of the real wave function qrit) for the Ising states r. The time dependent probability 
distribution of a generalized Ising model obtains as Prit) — q^{t). The functional integral employs 
a lattice regularization for single Weyl or Majorana spinors. We further introduce the complex 
structure characteristic for quantum mechanics. Probability distributions of the Ising model which 
correspond to one or many propagating fermions are discussed explicitly. Expectation values of 
observables can be computed equivalently in the classical statistical Ising model or in the quantum 
field theory for fermions. 



I. INTRODUCTION 



It has recently been shown that quantum mechanics can 
be obtained from a classical statistical ensemble [1]. This 
I— requires a suitable law for the time evolution of the proba- 
1^ bility distribution Pt (t) for the states r of the classical sta- 
I tistical ensemble, and an appropriate selection of observ- 
^*bles and their correlations describing sequences of mea- 
^ ^urements. An important criterion for the selection of the 
i__appropriate evolution law and observables is the compat- 
ibility with a "coarse graining" of the information to an 
isolated subsystem, for which only part of the information 



contained in {p^(t)} is available and used. 



00 The general setting how familiar "no go theorems" based 
•/^Qn Bell's inequalities 0; or the Kochen-Specker theo- 
rem [J] are circumvented has been discussed in ,1] both on 
^Nl fhe abstract level and by simple concrete examples. For a 
practical use of the equivalence between quantum mechan- 
ics and a suitable classical statistical ensemble one needs an 
^^implementation for realistic systems, as a quantum particle 
^ in a potential or quantum field theories. We have presented 
k>! ^ classical probability distribution in phase space which 
rNcan describe a quantum particle in a potential, including 
j^all characteristic "quantum features" as interference in a 
double slit situation, tunneling and the discreteness of sta- 
tionary states [5j]. For this purpose one has to postulate a 
specific "unitary time evolution" for the probability density 
in phase space which replaces the Liouville equation and is 
no longer compatible with the notion of particles following 
classical trajectories. While the choice of this evolution law 
can be motivated by the experimental success of quantum 
mechanics, one would also like to establish simple criteria 
for the fundamental time evolution law for probabilities in 
a classical statistical ensemble which do not rely anymore 
on the notion of individual trajectories. We believe that 
the answer to this problem is rooted on a more fundamen- 
tal level, namely quantum field theory for many particle 
systems, from which the notion of a single particle in a 
potential emerges as a special case. 

For the description of a quantum particle by a classical 
statistical ensemble in phase space the position and mo- 



mentum observables are "statistical observables" . The ex- 
pectation values of statistical observables are computable 
from the probability distribution, despite the fact that sta- 
tistical observables have no fixed value in the "classical 
states" (which correspond to points in phase space in this 
case). While statistical observables are familiar in classical 
statistical physics - an example is the entropy - they are 
usually not employed on a microscopic level. This raises 
the interesting question if the statistical observables on the 
level of one particle states can be obtained on a more fun- 
damental level from classical observables which take a fixed 
value in every state r. It seems natural that this fundamen- 
tal level corresponds to a quantum field theory. Indeed, we 
have demonstrated [6] that statistical observables and non- 
commutative products of observables can arise naturally 
from classical observables on a more fundamental level. 

A natural goal for the embedding of quantum mechanics 
into a fundamental description by "classical" probabilities 
seems to be the explicit construction of a classical statis- 
tical ensemble whose probability distribution describes a 
quantum field theory. As a first step, we have constructed 
the "classical" probability distribution for a quantum field 
theory of free massless fermions in two dimensions . This 
employs classical states r which correspond to the states 
of Ising spins on a one-dimensional lattice or chain. For 
every point x one introduces an occupation number n{x) 
which can take the values one or zero. (Equivalently we 
could use the Ising spins s{x) = 2n{x) — 1 which have val- 
ues ±1.) A classical state r is a sequence of occupation 
numbers r = {n(a;)}. The analogy to occupation numbers 
for fermions is very direct. Indeed, one has a map between 
states r and basis elements gr of a real Grassmann algebra 
[^] and one can associate the probability distribution pr (t) 
to a Grassmann functional integral [7}. 

In this paper we take a second step and construct the 
probability distribution which describes a quantum field 
theory of massless Weyl or Majorana fermions in four di- 
mensions. We will generalize this setting to Dirac spinors in 
an electromagnetic field in a forthcoming publication, but 
already the case of massless fermions has to circumvent a 
few obstacles and provides for important lessons. The ba- 
sic setting is a generalization of ref. 0] . The space points 
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X are now on a cubic three dimensional lattice, and we em- 
ploy fom- different occupation numbers n^(x),^ = 1 . . .4. 
The states of the classical ensemble are given by sequences 
of occupation numbers, t = {n^(x)}. 

An important ingredient for the formulation of a fun- 
damental evolution law for the probability distribution 
{Prii)} is the "classical" wave function qrit) d, For 
every time t this associates a real number g,- to every state 
r. The probabilities Pr are the squares of qr and therefore 
automatically positive, Pt = q^. The normalization of the 
probability distribution requires 

r T 

In turn, the wave function can be obtained from the prob- 
ability distribution up to a sign s,- = ±1, g,- = Sry^. For 
all quantities that can be computed from {pr} the choice 
of the signs corresponds to a choice of gauge without ob- 
servable consequences. (Convenient gauge choices respect 
continuity and differentiability of qr (t) , which excludes ar- 
bitrary "jumps" of St and fixes Sr to a large extent by the 
properties of {pr}- Then {qr} is essentially coinputable 
from {pt}, up to a few remaining "overall signs" [5|.) The 
concept of the classical wave function permits the formu- 
lation of simple time evolution laws Q which preserve the 
normalization of the probabilities. It is sufficient that the 
time evolution is described by a rotation of the real vector 

St, 

p 

The simplest evolution laws are linear - the rotation matrix 
Rrpit'^t) does not depend on the wave function {qrit)}. 

The classical wave function {qrit)} is a property of the 
classical statistical ensemble and resembles in certain as- 
pects to the Hilbert space formulation of classical mechan- 
ics by Koopman and von Neumann [9|. In contrast to ref. 
it is, however, a real function. This avoids the intro- 
duction of additional degrees of freedom which would cor- 
respond to the phases of a complex wave function. On the 
other hand, the physics of phases is one of the most char- 
acteristic features of quantum physics. We therefore have 
to implement a complex structure which maps the real 
classical wave function {qrit)} to an associated complex 
quantum wave function. For four dimensional fermions a 
natural complex structure can be associated to the equiva- 
lence between Majorana and Weyl spinors [U. Weyl 
spinors are directly formulated as complex entities, and 
the wave function for a single propagating Weyl fermion 
is indeed a complex function. On the other hand, a real 
representation of the Clifford algebra can be employed for 
Majorana fermions. The real wave function for a single 
particle can then be mapped to the complex wave function 
for Weyl spinors. 

We will extract the real wave function {qrit)} and its 
time evolution from a Grassmann functional integral based 
on a real Grassmann algebra. The action will involve four 
"real" Grassmann fields '07(^1 2;), but no separate conjugate 



fields ipj{t,x) or ■ipj{t,x). So far, the construction of the 
classical wave function from a Grassmann functional inte- 
gral has been based on pairs of conjugate Grassmann vari- 
ables ip and tp 0] ■ We have to generalize this construction 
for the case where no conjugate fields ip are available. Fur- 
thermore, a well defined setting requires a regularization of 
the functional integral that we implement on a space-time 
lattice. Since we want to be able to describe single prop- 
agating Weyl or Majorana fermions we have to develop 
a discretization for single Weyl spinors without "lattice 
doublers" . These constructions are all performed explic- 
itly and we end with a real Grassmann functional integral 
which describes a single species of massless Weyl or Majo- 
rana spinors in Minkowski space. We explicitly infer the 
wave functions for multi-particle and hole states. As usual 
in quantum mechanics, they depend on "initial conditions" 
set by the wave function {qrit in} a same initial time ti„. 

For a given action specifying the quantum field theory 
our construction allows the explicit computation of the 
quantum wave function {qrit)} for an arbitrary initial wave 
function {(/r(im)}- We may therefore ask what are the pos- 
sibilities for static wave functions that are invariant under 
Poincare transformations in the continuum limit. Any such 
wave function characterizes a possible "vacuum state" . We 
find by explicit construction that the vacuum state is not 
unique for a quantum field theory of free massless fermions. 
For every possible vacuum state we can construct one- 
particle and one-hole excitations corresponding to single 
propagating fermions. Their wave function obeys a free 
Dirac equation. Also the construction of multi-fermion 
states by applying appropriate annihilation and creation 
operators to the vacuum state is the same for all possible 
vacuum states. 

The physical meaning of the degeneracy of vacuum states 
has not yet been eludicated. It is an interesting speculation 
that they cannot be distinguished by "macroscopic observ- 
ables". In this case a "coarse graining" which averages 
over different vacuum states (and the associated excita- 
tions) would seem appropriate. Some light may be shed on 
this question by an investigation of massive Dirac spinors in 
an external electromagnetic field. Massless Dirac spinors 
can be understood as two species of massless Majorana 
spinors, but we do not know at the present stage the effect 
of the mass and the external potential on the degeneracy 
of the vacuum or more general ground states. 

For any given vacuum state we construct the one-particle 
wave function and discuss explicitly the corresponding 
probability distribution of the generalized Ising model. The 
time evolution of the wave function obeys the relativistic 
Dirac equation for a free massless Majorana particle or an 
equivalent Weyl-fermion. We specify observables for the 
position of the particle both as classical observables for the 
Ising-type model and as Grassmann operators. Expecta- 
tion values of those observables can be computed equiva- 
lently from the classical statistical ensemble for the gener- 
alized Ising model, or from the Grassmann functional of the 
quantum field theory for fermions. We also discuss general- 
ized one-particle states. For appropriate versions both mo- 
mentum and position of a single particle can be formulated 
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as classical observables. We introduce the non-commuting 
product of observables which underlies the measurement 
correlations of arbitrary functions of position and momen- 
tum. 

This paper is organized as follows. We introduce in sect. 
mi the action as an element of a real Grassmann algebra 
based on "real" Grassmann variables ipj{t,x). We employ 
a lattice discretization and establish the Lorentz-symmetry 
of the continuum limit. We also define a regularized Grass- 
mann functional integral. Sect. IIIII discusses the map be- 
tween this setting and the probability distribution of a gen- 
eralized Ising model. The Grassmann functional integral 
permits to compute the probability distribution {prit)} for 
every time t. It obtains by "integrating out" the Grass- 
mann variables for t' > t and t' <t. 

In sect. IIVI we turn to the time evolution of {prit)} and 
the associated "classical wave function" {qrit)} which is 
induced by the action formulated in sect. |lTl We show 
that it is "unitary" , corresponding to a rotation ^ of the 
classical wave function with R independent of {q}. We also 
investigate the notion of a conjugate wave function and use 
this to establish for the "partition function" the normaliza- 
tion Z = \. In the continuum limit the evolution equation 
becomes a type of Schrodinger equation for the real wave 
function {<Zt(0} which leaves the norm of the wave function 
invariant. On the level of the classical probabilities Prit) 
the evolution is non-linear. We also establish the evolution 
equation for the associated Grassmann wave function. In 
sect. |V]we discuss the classical observables of the general- 
ized Ising model in the usual setting of a classical statistical 
ensemble. They can be mapped to associated Grassmann 
operators, and the time evolution of expectation values can 
be described by a type of Heisenberg evolution equation. 
In particular, we discuss in this language the conserved 
quantities of our Ising type model. 

In sect. IVII we turn to the interpretation of our system in 
terms of particles. The Ising model describes an arbitrary 
number of massless free propagating Majorana spinors or 
the associated holes. One particle states obey a Lorentz 
covariant Dirac equation. We translate this description 
of Majorana spinors to an equivalent description of Weyl 
spinors. In this language the one-particle wave function be- 
comes complex. The general solutions for the one-particle 
or one-hole states are most easily discussed in this setting. 
The generalization to multi-fermion states is simple and 
leads to the usual totally antisymmetric wave functions for 
fermions. We discuss the associated symmetries, including 
continuous chiral symmetries and parity. 

In sect. IVIIII we discuss more systematically the com- 
plex structure which is associated to the equivalence be- 
tween Majorana and Weyl spinors. We group the four 
"real" Grassmann variables into two complex Grass- 
mann variables C,a ■ This maps the real Grassmann algebra 
constructed from ■0^ to an associated complex Grassmann 
algebra based on Cq. The transformation between -0^ and 
(Cai Ca) is a complex similarity transformation. Real Grass- 
mann operators built from and d/d^^ are mapped to 
complex Grassmann operators constructed from C,a^ d/d(^a 
etc.. Employing this complex structure, the action becomes 



the familiar complex action for free Weyl spinors. 

In sect. llXI we address the particle-hole conjugation. On 
the level of the classical Ising model this corresponds to 
an exchange of occupied and empty bits, or to a flip of 
sign of the Ising spins s-y(x). The particle-hole conjuga- 
tion is realized by a map within the Grassmann algebra 
which is closely related to the exchange of conjugate ele- 
ments. (It should not be associated to the usual charge 
conjugation in particle physics: Majorana spinors are in- 
variant under charge conjugation or simply flip sign.) The 
particle-hole conjugation maps states with n particles into 
states with n holes. In particular, we have investigated pos- 
sible vacuum states that are static, Lorentz-invariant and 
translation invariant and also symmetric under particle- 
hole conjugation. We show that with these conditions the 
vacuum state is not unique. In this section we also discuss 
position and momentum observables for various versions of 
one-particle states. Conclusions are presented in sect. [X] 



II. QUANTUM FIELD THEORY FOR FERMIONS 
IN FOUR DIMENSIONS 

In this section we formulate the quantum field theory 
for free Majorana spinors in four dimensions in terms of a 
Grassmann functional based on a real Grassmann algebra. 
The functional integral is regularized on a lattice. 

1. Action 

Let us consider the action 



(3) 



It involves four Grassmann functions ^-y(t,a:),7 = 
1 ... 4, a; = (xi, X2, X3). The integral extends over three di- 
mensional space and time, with dt = d/dt and dk — d/dxk- 
The real symmetric matrices Tk = (Tk)"^ are given by 



Ti = 



/O, 0, 1, 0' 

0, 0, 0, 1 

1, 0, 0, 

Vo, 1, 0, 0, 

1, 0, 0, 

0, 1, 0, 

0, 0, -1, 

0, 0, 0, -1 



/ 0, 0, 0, 1 

0, 0, -1, 

0, -1, 0, 

V 1, 0, 0, 



(4) 



and we deal with a real Grassmann algebra. Summation 
over repeated indices is implied. Within the Grassmann 
algebra the operation of transposition amounts to a total 
reordering of all Grassmann variables. The action ([3]) is 
antisymmetric under this operation. 



-S. 



(5) 



If we define formally the "Minkowski action" Sm = iS, the 
latter is hermitean, Sm — since S^j — —Sm- We will, 
however, not use the Minkowski action for our formulation. 
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2. Lorentz symmetry 

The action (|3]) is invariant with respect to four dimen- 
sional Lorentz-transformations, as defined infinitesimaUy 

by ip^it,x) -> 'ip'^{t,x) -\- S^p^{t,x), with e^^ = -e,.^,/i = 
(0,fc), 



(6) 



and ip!y{t,x) = ^j{t,x) — {^^dt + ^''dk)'(pj{t, x) accounting 
for the Lorentz transformed time and space coordinates. 
The matrices S^'' are given by 



^Ofe 



with 



0, -1, 0, 

1, 0, 0, 
0, 0, 0, -1 
0, 0, 1, 



1 



klra jrp 



(8) 



The Lorentz invariance of the action S is most easily 
established by employing the real matrices 



/ 0, 0, 0, 1 

0, 0, 1, 

0, -1, 0, 

V-1, 0, 0, 



(9) 



such that 



S 



(10) 



where 9o = dt and = i^s{l^)s'y (We consider if} here 
as a vector with components ipry and suppress the vector 
indices.) The real 4x4 Dirac matrices obey the Clifford 
algebra 



(11) 



with signature of the metric given by 77^^ = 
diag{— 1,1^1,1). This can be easily verified by using 
the relations 



{Tk,Ti} = 26ki , {7°,rfc} = 0. 
Furthermore, one finds 

(7T = -7V (7T=7^ 
and the relations 

[Tk,Ti] = 2ekimiT^ , [TkJ] =0 , 
7V7V = I, {7",n-0, {7',/}-0. 



-1, 



(12) 
(13) 

(14) 



The Lorentz generators ([7|) obtain from the Dirac matrices 
as 



(15) 



With 

one finds that ^ transforms as 



(16) 



(17) 



such that ipTp = '07'07 is & Lorentz scalar. We recognize 
in eq. (|10l) the standard Lorentz invariant action for free 
fermions in a Majorana representation of the Clifford alge- 
bra with real 7^-matrices. 



(7) 3. Discretization and functional integral 

The action ^ may be considered as the continuum limit 
of the regularized action 



with 



5 = E 



L{t) = ^il}j{t,x)B^{t + e,x). 



(18) 



(19) 



We will place x on points of a cubic lattice and consider a 
regularization where e is small compared to the lattice dis- 
tance, permitting a continuum limit e -> for time even for 
discrete space points. Here B^{t, x) is a linear combination 
of Grassmann variables iplt^x), which we write as 

B^{t,x) = 'ip^{t,x) - eF^[ip{t);x], (20) 

such that 



The sum in eq. p8|) extends over discrete time points t„, 
with = eJ2t = ^J2n ' - in = e,n € Un < i„ < 
tf. We employ 



9i^(t)--[^(t + e)-V(i)]: 



(22) 



such that the Grassmann property ipy{t, a;) = results in 

'ipj{t,x)dti'y{t,x) = - ^j{t, x)^jjj{t + e, x). (23) 

The continuum limit is taken as e — >■ for fixed tin,tf. 

Similarly, we sum in eq. (PT|) over points cc of a cubic 
lattice with lattice distance A and = ^^J2x- Corre- 
spondingly, is defined as 

[^{t + e); x] = {n)^sdkMt + e,x)+ 0(e) (24) 

with use of the lattice derivative 

dk^dk-\j2 ^ki^TiTj-iHk . (25) 
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Here we employ 

4^7 W ^ ^[^^{x + Ak) - iJ^{x ^ Ak)] (26) 

and 

hM^) = ^{tp-,{x + Ak)+ijj{x^Ak)-2^^{x)}, (27) 

with Afc the unit lattice vector of length A in the k- 
direction, A| — A^ . Details of the choice of the discretiza- 
tion and properties of the lattice derivative are discussed in 
appendix A. For our purpose it is important that the dis- 
cretization does not introduce any fermion doubling, such 
that we can describe single Weyl or Majorana spinors. The 
term ~ 6k, which prevents the occurrence of doublers, van- 
ishes in the continuum limit A — )■ 0. 

In appendix A we also specify the omitted term ^ 0(e) in 
eq. (|24p which vanishes for e -> or A 0. This term is 
needed in order to guarantee that the Grassmann element 
By(t,x) has the same normalization as ipy{t,x). In other 
words, we can obtain from by a rotation 

Bfi^) = ^Rfs{x,y)ips{y), 
y 

R-rjsiz, y)R^s{x, y) = 5n-,5{z, x). (28) 

V 

We may take the limit e — first. At this stage the model 
is formulated on a discrete space lattice. The continuum 
limit in the space direction, A — >■ 0, can be taken at the 
end. We observe that the space-continuum limit, leading to 
the action ([3]), also involves a rescaling of the Grassmann 
variables by a factor A'^"^/^-'. The Lorentz symmetry is re- 
alized in the continuum limit where both e and A approach 
zero. 

The functional integral is defined by the partition func- 
tion 

Z = j V^gf [^{tf)] e-^g,n [^{Un)] , (29) 
with the functional measure 

f 1^^ = 11 [ {dMt,x)...dMt,x)). (30) 

The boundary terms and gf only depend on the Grass- 
mann variables V'(ijn) ^^^id 'ip{tf ), respectively. As we will 
see below, the boundary terms gf and gin are related to 
each other, such that the functional integral (I^Hl) is fully 
specified by the choice of gin ■ 



III. PROBABILITY DISTRIBUTION AND WAVE 
FUNCTION 

In this section we associate to the Grassmann functional 
integral (|29p a family of classical probability distributions 
{Pr{i)}- For every given time t this probability distribution 



associates to each classical state r a classical probability 
Prit), which is positive and normalized 

Pr(t)>0, Y.pr{t) = l. (31) 

r 

The positivity and normalization of the probabilities holds 
for an arbitrary choice of g™, provided ^/ is related to gin 
appropriately. For every given gin the probability distribu- 
tion {prit)} is uniquely computable for all t, such that the 
functional integral also specifies the time evolution of 
the probability distribution. 

The states t of the classical statistical ensemble are dis- 
crete if we use the formulation of the functional integral on 
a space lattice. If we assume a cubic lattice with lat- 
tice points we associate to every point x four bits (or Ising 
spins) nj{x) that can take the values 1 or 0. The states r 
are then given by the sequences of bits, r = {n^(a;)}, each 
T corresponding to an ordered chain of AL^ numbers 1 or 

0. There are a total of 2**^ different states t. 

1. Grassmann wave function 

The construction of the probability distribution {pr} will 
rely on the isomorphism between states r and the basis el- 
ements gr of a Grassmann algebra that can be constructed 
from the Grassmann variables ip-yix) for fixed t. Each basis 
element gr is a product of Grassmann variables 

gr ^ ip.y,{xi)ipy^{x2) ■ . ■ (32) 

which is ordered in some convenient way. To be specific, we 
define some linear ordering of the lattice points and place 
variables with "smaller x" to the left, and for each x place 
smaller 7 to the left. If a Grassmann element gr contains a 
given variable ip-y (x) we may put the number (x) in the 
sequence t to 0, while we take n^{x) = 1 if the variable 
ip-f{x) does not appear in the product p2p . This specifies 
the map between the 2"^^^ independent basis elements gr 
of the Grassmann algebra and the states r. 

Our first step computes for the functional integral ([^ a 
Grassmann wave function g{t), which is an element of the 
Grassmann algebra constructed from the Grassmann vari- 
ables "07(^7 x) at given t. In order to define this Grassmann 
wave function we decompose the action (|18p 

S = 5<+5>, 

5< = 5]L(0, Sy=Y^L{t'). (33) 

t'<t t'>t 

The wave function g(t) obtains by integrating out all Grass- 
mann variables for t' < t 

git) = J ViPit' < t)e-^<g,n- (34) 

We observe that g{t) depends on the Grassmann variables 
^{t). More precisely, it can be constructed from linear com- 
binations of products of the Grassmann variables ^^{t,x) 
for fixed t. 
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2. Classical probabilities and wave function 

We next expand g in terms of the basis elements gr of 
this algebra 

5W -E'?-(*)-94V'W]- (35) 

r 

We would like to associate the real coefficients qr{t) with 
the classical wave function, such that the classical prob- 
abilities obtain as Prii) = Iri^)- This requires for ev- 
ery t the normalization X]t9t(^) ~ 1- We will show in 
the next section that this normalization condition is in- 
deed obeyed, provided it holds for the initial wave function 

The conjugate wave function is defined as 

J V^{t' >t)gfe~^>. (36) 

Again, this is an element of the Grassmann algebra con- 
structed from V'(i)- In terms of g and g the partition func- 
tion reads 

j vmmgit). (37) 

where the Grassmann integration J 'D'ijj{t) extends only 
over the Grassmann variables ^^{t^x) for a given time t. 
The conjugate basis elements of the Grassmann algebra gr 
are defined by the relation 

~9r9r ^XlXl^^ix) (38) 

X 7 

(no sum over r) and the requirement that no variable ^^7(2^) 
appears both in gr and gr- They obey 

j vm~9r [m] 9p [m] = Srp. m 

Expanding 

5(t)=^g,(i)5r[V'(i)] (40) 

r 

yields 

Z = Y^qr{t)qr{t) (41) 

T 

We will see in the next section that for a suitable choice of 
gf the conjugate wave function obeys for all t the relation 
qr{t) — qr{t). Together with the normalization condition 
J^r'lr ~ 1 ^^^^ guarantees Z = 1. 

In consequence, we can express the classical probabilities 
Prit) directly in terms of the Grassmann functional integral 

Pr{t) = J Vi:gfVr{t)e'^g,n, (42) 

with Vr{t) a projection operator 

Vr{t)gp[iAt)\ = gr[^{t)\5rp. 

J V^P{tyg4^j{t)]Vr{t)gp[m] ^ Sra5rp (43) 



The formal expression of Vr in terms of the Grassmann 
variables V'(^) and derivatives d/d'ip{t) can be found in ref. 
[6]. 

The interpretation of the Grassmann functional inte- 
gral in terms of classical probabilities is based on the 
map g{t) — > {^^-(t)}, which in turn is related to the map 
{Qt} {Pt} ~ {<lr}- The opposite direction of a map 
{pr{t)} — > g{t) is also possible, provided we choose an ap- 
propriate sign convention for s^- = ±1 in 

qr{t) = Sr{t),/Mt)- (44) 

This sign convention corresponds to a choice of gauge 
which does not matter if all expectation values and cor- 
relations can finally be expressed in terms of the proba- 
bilities {pr{t)}. The map g{t) ^ {qrit)} is invertible. A 
given "quantum state" or equivalent classical statistical en- 
semble may be specified by the "initial value" at some time 
to,g{to), or the associated wave function {qr{to)} or classi- 
cal probability distribution {prito)}- This is equivalent to 
the specification of gm = g{tm) in the functional integral. 

IV. TIME EVOLUTION 

In this section we compute the time evolution of the 
wave function {qr{t)} and the associated probability dis- 
tribution {pr{t)} — {qri^)}- This will lead to a type of 
Schrodinger equation for the real wave function {qr{t)^, 
as well as an associated evolution equation for the Grass- 
mann wave function g{t). We will use the properties of 
this evolution equation in order to establish that the norm 
J2r 9r(0 is conserved, such that {g^(t)} can indeed be in- 
terpreted as a time dependent probability distribution. 

Due to the particular form of the action ^ , which only 
involves one type of Grassmann variables (and no conju- 
gate variables as in ref. [7] ) the discrete formulation of the 
functional integral is crucial. We will see that g{t) jumps 
between two neighboring time points, while it is smooth 
between g{t + 2e) and g{t). We will therefore distinguish 
between even and odd time points and use the definition 
([55)) of the wave function only for even times. For odd times 
we employ a different definition, which will guarantee the 
smoothness of the time evolution of {qrit)} for both even 
and odd time points. In the continuum limit, e — >■ 0, the 
time evolution of the wave function is then described by a 
continuous rotation of the vector {qr{t)}. In Appendix B 
we show that by a distinction between even and odd time 
points, with i}}^{t,x) = ip^it + e,x) for t even, we can ac- 
tually map our formalism to the formalism of ref. [7] with 
conjugate spinors ip. 

The time evolution of the Grassmann wave function 
obeys 

g{t + e)= Jvme-'^^'^git). (45) 

This determines {qr{t + e)} in terms of {qr{t)}. Thus the 
action (PT|) specifies the dynamics how the proba- 

bility distribution {pr{t)} evolves in time. The particular 
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dynamics of a given model is determined by the term ~ eF-y 
in eq. (piT) . This will play the role of the Hamiltonian in 
quantum mechanics. 

1. Trivial time evolution 

In order to understand the structure of the evolution 
law (|45l) we first consider the case F-y = 0. For a vanishing 
Hamiltonian this should describe a trivial time evolution 
with a static probability distribution. For F^ — the evo- 
lution equation (|45l) reads 

g{t + e) = J PVWexp{-^^^^(i,x)^^(i + e,x)}5(t). 

(46) 

We may write g{t) in a product form 

= n n ^) + a:)] , (47) 

with some fixed ordering convention of the factors assumed, 
e.g. smaller 7 to the left for given x, and some linear 
ordering of the lattice points, with "lower" points to the 
left. In the product form eq. yields 

g{t + e)^ [vmUIi (48) 

I [1 - tp-yit, x)^p-y(t + e, a:)] [aj{t, x) + b-y{t, x)'ip-y{t, x)] | 
= 1^7(^,2;) +r]yaj{t,x)->p^{t + e,x)y 

X 7 

Here we use the fact that each individual Grassmann inte- 
gration J d4!-f{t,x) can be performed easily. 



dip{l — 'ipip){a + bip) = b — aip, 



(49) 



and = ±1 results from the anticommuting properties of 
the Grassmann variables (p, with 

'/I = 7?3 = 1 , ??2 = 7?4 = -1- (50) 
As a result, we can write 

9{t + €) = Y,<lr{t)Cgr[i'y{t + e,x)\, (51) 

r 

where Cqt obtains from gT by the following replacements: 

(i) for every factor ili^{x) in gr one has a factor 1 in Cgr] 

(ii) for every pair [x, 7) for which no ip-f{x) is present in gr 
one inserts a factor riy'ip^{x) in Cgr- The ordering of the 
factors rj~f'ip^{x) is the same as the ordering assumed in the 
product (j47p . This implies that we can indeed write the 
action of C on the product ([Tf|) as 

Cg[t) - n n 2:) + V^a^{t, x)^^{t, x)) . (52) 

X 7 

The conjugation operator C maps every basis element g^- 
into its conjugate element g-r up to a sign Ut = ±1, 



(53) 



Applying C twice on the product (|47|) multiplies each factor 
by 77-y. The factors r]j drop out due to the even number of 
minus signs, such that C is an involution, = 1, or 



9t = gr- 



(54) 



This has the consequence that g{t + 2e) obtains from g{t) 
by simply replacing the arguments "07 (i, x) by ■0^(i-|-2e, x), 
or 



g^(t + 2e) = g^(i). 



(55) 



The trivial static evolution of the wave function is realized 
between t and t + 2e. 

The jump between g and Cg for neighboring time points, 
g{t -I- e) = Cg{t), suggests the use of eq. ((35)) for the def- 
inition of the wave function |(7^(t)} only for "even time 
points" , namely those that obey t„ = <i„ + 2me, m S N. 
For odd time points, i„ — Un + (2m — l)e, we employ a 
different definition, namely 



g{t)^J2lrit)C9r[m]- 



(56) 



With this definition one finds that for F.y = the wave 
function {qrit)} is indeed time independent, Qrit) = 

Similarly, we use the definition pO)) for qrit) only for 
even t, while for odd t it is replaced by 



5(<) =^'7r(t)Cg40(i)]. 



(57) 



The map C acts similarly as C, with rj^ replaced by ij^ = 
—Tj-y, fji — fis — — l,'ij2 — — 1, and = 1. One observes 



Cgr = VrO-rgr, 



(58) 



with rjr — —1 ii gr contains an odd number of Grassmann 
variables, while rjr ^ 1 for an even number. Relation 
follows from = 1 or 



Cgr — CTrgr, 



(59) 



together with the property that C obtains from C by an 
additional minus sign for every ip-iactoT 



Cgr = VrCgr- 

The relation ([55)) implies 

V^CgrCgp = 5rp 



(60) 



(61) 



such that the identity ((IT)) holds both for t even and t odd. 
Employing the identity 



~g{t~e)^ / 2?V(i)3(i)e-^^*"^\ 



(62) 



one finds for F^ — 



9{t - = X! ^r{t)Cgr [ip.yit - e, x)] , (63) 
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and therefore a time independent qr{t). With g{tf) = g/ 
we can choose gf = [^(^/)] ^'-'^ ^ven, or 

cjf = ['^(*/)] Then the relation 

9t(0 = ^T{t) holds for all t. We conclude that for = 
the dynamics is trivial, with time independent q-r {t) = q-r (i) 
and therefore a time independent probability distribution 

{Pr}- 

2. Unitary time evolution 

We next turn to the dynamics which is described by the 
Lorentz-invariant action or the corresponding regular- 
ized form IH]), dU), with given by eq. In con- 
sequence, we have to replace in eq. (|46p ip^{t + e,x) by 
B^{t -\- e, x) according to eq. (PH)) 

+ (64) 

a; 7 

This implies (for even t) 

ait + - II 9r(t)C5r [B^{t + e, x)] . (65) 

T 

Since B^{x) is related to ipj{x) by a rotation (pS)) . it is 
straightforward to show that (7t-[^7(2;)] is also connected 
to gr ['07(2;)] by a rotation among the basis elements 

gr[B^{x)] = ^gp[ipj{x)]Rpr, 
p 

'^RrpRap = Sra- (66) 
P 

One infers 

git + e) = Y,<lrit)C9p[i^jit + e)]Rpr 

= J2lrit + ^)Cgr[i'^it + e)] (67) 

r 

with a rotated wave function 

qr{t + e) = Y,Rrpqp{t). (68) 

p 

An analogous argument leads to the same evolution (|68l) 
for t odd. 

Rotations preserve the length of the vector {qr} such 
that ^^qr{t) is independent of t. Choosing g^ = 
Y.Tlr{'tin)gT['4'{tin)] with J^t Iritin) = 1 One infers 
Y^rlrit) = 1 for ah t. Therefore {pr{t)} = {q^t)} has 
indeed for all t the properties of a probability distribu- 
tion, namely positivity of all Pr and the normalization 
J2tPt = ^- In analogy to quantum mechanics we call a 
time evolution which preserves the norm of {qrit)} a "uni- 
tary time evolution" . A unitary time evolution is crucial 
for the probability interpretation of the functional integral 
(Ei). 



3. Evolution of conjugate wave function 

We next want to show the relation (for t even) 

ff(t-e) =^g,(t)Ci?,p5p[0^(t-e,x)]. (69) 

T,P 

The definition (|57p of the conjugate wave function g, 

ff(^-e) =5l9r(i-e)Cg.[V'7(i-e,a:)], (70) 

r 

then imphes 

qr{t-e)=J2qp{t)Rpr, (71) 
p 

such that 

g,(t) =^i?,pqp(t-e). (72) 

p 

Comparing with eq. ()68p one infers that qrit) and qr(i) 
obey the same evolution equation. (By an analogous argu- 
ment eq. ([7^ also holds for t odd.) If q and q are equal for 
some particular time to, they will remain equal for all t. 
In order to show eq. ([M)) we employ eq. (p^ . 

5(i - e) = / T^mgit) exp { - ^ V'7(^ - a;)B^(t, a;)} 

= j 'Dtl){t)g{t) exp{ -^■0^(i-e,2;)^7'5(2^'y)V'5(^, y)} 
x,y 

= /"l?V(05(i)exp{-^B^(i-e,a:)V'7(i,a;)}, (73) 

a; 

which replaces in eq. ([M]) ■(/'7(^ ~ 2;) by i?^(t — e, x), 

B^{t,x) = ^■il;s{t,y)Rsj{y,x), 
y 

= Y.{R-%s{x,y)Mt,y). (74) 

V 

leading to 

g{t -t)=Y. 9(^)C3r [Mt - e, x)] . (75) 

T 

Eq. ([M)) then follows from the relation 

gr[B',{t,x)] ^^Rrpgp[4;j{t,x)]. (76) 
p 

Indeed, we can expand gT[B^{t,x)\ in terms of the com- 
plete set of basis elements cjp [i^^it, x)] , 

gr[B^y{t,x)] ^^Arpgp[ip^{t,x)], (77) 
p 
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with 

Arp - j V^tygr[B^{t,x)\gp[ii;^{t,x)\ 

= j Vi,{tyg,[^^{t,x)\gp[B^lt,x)\ (78) 

= ^ I 'Dtp{t)gr[^l)^{t,x)\ga[^^{t,x)\Rap ^ Rrp- 

The second hne uses the fact that RB = ip and the invari- 
ance of the functional measure / 'Dipit) under rotations 
ip — >■ Rip. This concludes the proof of eq. (p^ . 

4. Partition function 

If |(7r(to)} equals {^,-(^0)} for some time to we can use 
Qrit) = Qrit) for all t and infer from cq. (|^T|) 

^ = E9'w- (79) 

r 

As it should be, Z remains invariant under rotations (j68p of 
the vector {qt} and is therefore independent oft. Inversely, 
we may actually use the fact that the relation (|4T|) holds 
for all t in order to show that {qt} and {Qt} must obey the 
same time evolution provided q{t + e) obtains from q{t) by 
a rotation. Assume 

qrit + e)=J2^rpqpit) , (t + e) = ^ i.p^p (t ) (80) 
p p 

for some arbitrary regular matrices A and A. The time 
independence of eq. (j4ip then implies in a matrix notation 

Z = g^(t + e)q{t + e) = g^(i)i^Ag(i) = f{t)q{t). (81) 

This holds for arbitrary q and (7 (corresponding to arbitrary 
9in,gf) provided 

A^A=1 , A={A^)-\ (82) 

If ^ describes a rotation the matrices A and A coincide. 
The relation (|82p demonstrates that the time evolution of 

by a rotation (|68p is a sufficient and necessary condition 
for q{t) — q{t) and therefore for the expression p2)) of the 
probability in terms of the Grassmann functional integral. 
li A is not a rotation matrix, the relation q{t) — q{t) cannot 
be maintained for arbitrary t and arbitrary g(to). 

These remarks allow for straightforward generalizations 
of our setting. Whenever we can write 

L{t) ^^'ip-^{t,x)B^{t + e,x), (83) 

X 

with _B-y(t + e, a;) containing terms with an arbitrary odd 
number of Grassmann variables ip-y{t + e,x), then eq. ()66p 
is sufficient to guarantee q^ — const, and q{t) = q(t) 
provided g(to) = ^(io)- 



5. Boundary terms 

The final point we have to settle in order to establish the 
normalization Z = 1 and the expression (j42l) concerns the 
equality of 9(^0) smd q{to) for some arbitrary time ^o- This 
is achieved by a proper choice of the relation between the 
boundary terms gf and g^ in eq. (j29p . For this purpose we 
may imagine that we (formally) solve the evolution equa- 
tion (|68l) in order to compute { (^ / ) } in terms of { qr {tin } , 
9in = J2t QT{tin)gT \}p{tin)\ ■ Let US assumc that </ is even 
such that 

9{tf) = (84) 

r 

~9{tf) = Y.^r{tfygr[Htf)\ =9f- (85) 

T 

It is then sufhcicnt to choose gf such that qritf) — qritf). 
Equivalently, we may specify the wave function 1 5^.(^0)} = 
{9t(^o)} at some arbitrary time to and compute the cor- 
responding gin and by a solution of the evolution equa- 
tion, using the fact that the rotation (|68|) can be inverted 
in order to compute q(t — e) form q{t). 

6. Continuous evolution equation 

Finally, we cast the evolution law into the form of 
a differential time evolution equation by taking the limit 
e — > 0. We employ eq. (|65|) consecutively for two time 
steps 

dt9it) = ^[<?(i + 2e)-5(t)] 

= Y^^qT{t)\^gr[ilJ^{x) ~ 2eF.^{x)] -gr[ipt{x)]^ 

T 

= ^Qr{t)gp[!p^{x)]Kpr 

= J2dtqr{t)gr[i'^{x)], (86) 

r 

resulting in a Schrodinger type equation for the real wave 
function {qr(i)}, 

dtqrit)=J2K^pqpit). (87) 
p 

(In eq. ((86|) we use a fixed basis, corresponding to the 
basis elements gT constructed from •ip{t) for g{t), and from 
'ip{t + 2e) for g{t + 2e).) Since the evolution describes a 
rotation, the matrix K is antisymmetric 

Kpr = -Krp. (88) 

We identify this evolution equation with the Schrodinger 
equation for a quantum wave function for the special case 
of a real wave function and purely imaginary and hermitean 
Hamiltonian H — ihK. 

The time evolution (l87l) translates directly to the prob- 
abilities (no summation over r here) 

dtPr 2 E KrpSrSpy/prPp- (89) 

p 
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Once the signs Sr(io) are fixed by some appropriate con- 
vention at a given time to, the signs ST-it) are computable 
in terms of the probabihties. This follows since for all t the 
wave function Qrit) is uniquely fixed by qrito) or pT-ito), 
and Pr^t) is uniquely determined by qrit). In principle, it 
is therefore possible to formulate the time evolution law for 
the probabilities uniquely in terms of the probabilities. In 
particular, the positive roots ^fp^ obey 



dt\/P^ = ^ SrSpKrp^/p^. 



(90) 



This is a standard linear equation for ^/p^, except for the 
additional sign information stored in {sr(i)}- The role 
of the sign is to prevent yfp^ to become negative during 
the evolution and to maintain the normalization condition 
'Y^^Pt = 1- Whenever for a fixed sign Sr the evolution 
([ro]) would lead to crossing zero, the sign Sr flips 

such that yfp^ increases in the following times instead of 
decreasing to negative values. In principle, this provides 
for an algorithm of keeping track of the appropriate signs 
{sT-{ty\. Whenever Pr{t{)) = for some time t^, only one 
choice of ar {to + e) is compatible with p^ (to + e) > and 
'Y^Prito + e) = 1- This condition decides if CTt- jumps at 
to or not. 

The elegant way of keeping track of the sign informa- 
tion is, of course, the use of the wave function. The linear 
evolution law (|87l) is very convenient for the description 
of the classical ensembles which obey an evolution law of 
the type (lUUl) . More generally, the general form of the evo- 
lution law (|87p allows for the construction of simple time 
evolution laws for classical probability distributions. For 
a general evolution law the basic property to be respected 
is the condition of unit norm of the probability distribu- 
tion. This can be quite cumbersome for a general evolution 
equation for {pr}, but it is extremely simple on the level 
of {qr} where only the length of a real vector has to be 
preserved. 

7. Grassmann evolution equation 

The matrix K can be extracted from the Grassmann 
evolution equation 



dtg = ICg, 
with Grassmann operator /C given by 



d 



■{Tk)^sdkil^s{x). 



(91) 



(92) 



In order to proof the relation ([M)) we consider the 

Grassmann operator 



d d 



d'lp-y{x) 

1 - 2eF^{x) 



d 



d^j{x) 



(93) 



Its action on an arbitrary Grassmann element g results 
in the replacement of the particular Grassmann variable 



■0-y(a;) by ip-yi^) ~ 2ei^^(x). In terms of these operators we 
can write 

git + 2e) = {J[ W^{x))git) ^ U{t + 2e, t)g(t), (94) 



with U the Grassmann evolution operator [3] , [6| . (In eq. 
(|94p we only have kept contributions up to linear order in 
e.) We can now verify the relations 



^/C^(a 



x,7 



d 



W^{x) = exp{2e/C^(x)}, 
such that eq. (|94|) results in 

g{t + 2e) = exp(2e/C)5(t) = g{t) + 2e]Cg{t), 



(95) 



(96) 



thus establishing eq. ([?T|) . (Again, at several steps we have 
neglected corrections ^ that vanish in the continuum 
limit.) 

From the expansion of K,gr in the basis elements gp, 



l^gr ^^gpKpr, 



we infer the matrix element 

= / T>ijjgplCgr 



(97) 



(98) 



Taking finally the space-continuum limit by rescaling ip-y (x) 
and d/d^j{x) such that 



d 



dlPryi^X) 



i'siy) } = SjsS^ix - y), 



(99) 



we arrive at the continuum form of the Grassmann evolu- 
tion equation 



dtg ^ICg , /C = 



d 



X dipj{x 



■{Tk)jsdkil^5ix) 



(100) 



This evolution equation will be the basis for the interpre- 
tation of the time dependent wave function {qr(i)} and 
probability distribution {pr jt)} in terms of propagating 
fermionic particles in section IVll 



V. OBSERVABLES 

Classical observables A take a fixed value Ar for every 
classical state r. In classical statistics the possible out- 
comes of measurements of A correspond to the spectrum 
of possible values Ar . The expectation value of A obeys 



(A) ^Y^PrAr. 



(101) 
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Our description of the system will be based on these clas- 
sical statistical rules. For example, we may consider the 
observable measuring the occupation number N^{x) of the 
bit 7 located at x. The spectrum of possible outcomes of 
measurements consists of values 1 or 0, depending if a given 
state T = has this particular bit occupied or empty. 

1. Grassmann operators 

For the Grassmann basis element gr associated to t 
one finds Nj{x) = if g^- contains a factor ij;j{x), and 
Nj{x) = 1 otherwise. We can associate to this observable 
a Grassmann operator Af-y (x) obeying (no summation here) 



Afj{x)gr = {Ny{x))^gr , JV^{x) = 



d 



dip^{x) 



ij^ix). (102) 



Two occupation number operators Af-yj^{xi) and Af-y2{x2) 
commute if 71 7^ 72 or xi ^ X2- 

In general, we may associate to each classical observ- 
able A a diagonal quantum operator A acting on the wave 
function, defined by 

{Aq)r = Arqr. (103) 

This yields the quantum rule for expectation values 

{A) = {qAq) =Y,qrArpqp, (104) 



with A a diagonal operator Arp ~ ArSrp- In the Grass- 
mann formulation one uses the associated Grassmann op- 
erator A obeying 



Agr = Argr 



such that 



V^gAg. 



(105) 



(106) 



Here g is conjugate to g, i.e. for g — J2t 1^9^ o^i^ has 

2. Time evolution of expectation values 

In classical statistics the time evolution of the expecta- 
tion value is induced by the time evolution of the proba- 
bility distribution 



(107) 



This corresponds to the Schrodinger picture in quantum 
mechanics 



{A{t)) = {q{t)Aq{t)), 



(108) 



or the corresponding expression in terms of the Grassmann 
algebra 



(Ait)) = / Vil;g{t)Mt)- 



(109) 



Using dtq = Kq (|87p we infer for the time evolution of 
the expectation value the standard relation 



dM) = {q[A,K]q). 

Similarly, we use dtg ^ tCg ([9T|) and 
dtg = -IC^g, 



(110) 



(111) 



with K.^ obeying for arbitrary Grassmann elements / and 
g the relation 



Vi^lC^-gf = I Vr9lCf. 



One infers 



dt{A) = Jvi;g[A,JC]g. 



(112) 



(113) 



3. Conserved quantities 

Conserved quantities correspond to Grassmann opera- 
tors A that commute with JC. Let us consider operators of 
the type 



d 



dtpeiy) 



^er,(2/)'0,,(y), 



(114) 



with B^ri depending on y and derivatives with respect to y. 
They obey the commutator relation 



(115) 



d 

where all quantities on the r.h.s. depend on yk and dk ~ 
d/dyk- We can use eq. (|115p in order to find the conserved 
quantities 

Bi = I -^6i(9fc)V7 , ^2 = / ■^i^sh2{dk)i^s. (116) 

Here we use the property that / commutes with Tk and 
&i_2 contain only derivatives with respect to y. Hermitean 
(symmetric) operators are found if bi only involves even 
numbers of derivatives and 62 only odd numbers. We ob- 
serve the presence of off-diagonal conserved quantities as 



(117) 



which are not classical statistical observables with fixed 
values in every state t, but can nevertheless be useful for 
the understanding of the system. 
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VI. PARTICLE STATES 

Our system admits a conserved particle nmxiber, corre- 
sponding to the Grassmann operator J\f, 



d 



y dtpjiy) 



^^(y) , [AA,/C] -0. (118) 



The particle number is Lorentz invariant. We can decom- 
pose an arbitrary Grassmann element into eigenstates of 



— ^ ^ Am dm J A/" g-m — 



mg„ 



(119) 



The time evolution does not mix sectors with different par- 
ticle number m, such that the coefhcients Am are time in- 
dependent 

dtg = ^ Amdtgni , dtgm = K^gm- (120) 

m 

We can restrict our discussion to eigenstates of N . 
1. Vacuum 

Let us consider some static vacuum state go with a fixed 
particle number toq, 



/C.go = , TVgo = m^go , / 'Dikg^go = 1 



(121) 



An example for a possible vacuum state is the totally empty 
state 50 — |0), with 

io) = = n 11^^7(2^) - n('^i^2V'3V'4), 



X 7 



= 0. 



It obeys 



Vil)\Q) = 1 



1. 



(122) 



(123) 



Another example is the totally occupied state g^) — 1, with 
Too = B. We shift the particle number by an additive 
"renormalization" n = m — toq ? such that the vacuum cor- 
responds to n = 0, and g = A„(;„. An eigenstate of M with 
eigenvalue m = mo -(- n is called a n-particle state if n is 
positive, and a n-hole state for negative n. 

We recall that in a formulation with discrete time steps 
g{t) refers to even t. For odd t and continuous time evo- 
lution the Grassmann wave function is given by g{todd) ~ 
Cg{t). The conjugation C interchanges the role of particles 
and holes. Thus for g{t) — |0) one has g{todd) = |1) = 1, 
such that the vacuum Grassmann wave function flips be- 
tween fully empty and fully occupied for t even and odd. 
Similarly, the excitations flip between particles and holes. 
Possible vacuum states with the same Grassmann wave 
function for t even and odd, as go — (|0) + |l))\/2, will 
be discussed in sect. IIXI 



2. One-particle and one-hole states 



We next define creation and and annihilation operators 
d 



at(a;), a^{x) as 



al{x)g 



9 , a^{x)g = tp-i{x)g. 



(124) 



They obey the standard (anti-)commutation relations 
{a\^{x), ae(j/)} = 5^c5{x - y) , TV = / a^^{x)a^{x), 



[a\{x),N] -a\{x) , [a^{x),N] = a^{x). (125) 

Acting with the creation operator on the vacuum produces 
one-particle states 

9i{t) = / (l'y{t,x)a\{x)go = Gigo, 

J X 

{M~mo)gi=gi, (126) 

while a one- hole state with n = — 1 obtains by employing 
the annihilation operator 



9-iit) 



qy{t,x)a^{x)go 

qy{t,x){j°)jsas{x)go = G^igo, 
- mo)g-i ^ ~g-i, (127) 



with 



(128) 



If we transform the one-particle wave function q^(t,x) 
infinitesimally according to 



5q^ 



(129) 



the operator Qi is Lorentz invariant. (We omit here the 
part resulting from the change of coordinates.) Similarly, 
q^ transforms as q-y and the corresponding infinitesimal 
transformation 



(130) 



results in an invariant G-i- If the vacuum go is Lorentz 
invariant, the Lorentz transformations (|129p . (|130p of the 
one-particle or one-hole wave functions will obey the same 
evolution equations as the original wave functions. 

The time evolution of the one particle wave function q^ 
is given by 



d,g,^JCg,^ j{d,q^)go^ j^. 

and similar for the hole. With 

d 









J X 


'ic,— 


50 



d'ljj-y{x) 



£7 



(131) 



(132) 
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and 



(133) 



one obtains Dirac equations for real wave functions (with 
do = dt) 



rd^q = , rd^^q = , irfd^q = 0. 



(134) 



We emphasize that these equations foUow for arbitrary 
static states go which obey ICgo = 0. 

If needed, we may multiply Qi with an appropriate nor- 
malization factor such that the wave function obeys 



(135) 



and gi has a standard normalization. (No such factor is 
needed for go — |0).) We observe that the existence of 
particles and / or holes depends on the vacuum state go- 
For the example 5o = |0) no hole states are present, while 
for go = 1 no particle states are allowed. 

3. Weyl and Majorana spinors 

Let us consider go = |0) where a-y(x)|0) = implies that 
no hole states exist. There are then only particle states 
and the propagating degrees of freedom correspond to Ma- 
jorana fermions. In four dimensions Majorana spinors are 
equivalent to Weyl spinors flo'|. Indeed, we may introduce 
a complex structure by defining a two-component complex 
spinor 



(ySi = (7i + iq2 , </32 = 93 + ^qi■ (136) 



The matrices — 7^7'^ are compatible with this com- 
plex structure. They are translated to the complex Pauli 
matrices 



(137) 



while the operation (/?—>■ corresponds to the matrix 
multiplication q — > Iq. The operator Vk in eq. (|117p acts 
on the one-particle states as the usual momentum operator 

Vk"^ = ~idk(p. (138) 

The Dirac equation reads in the complex basis 

dtip = Tkdk(p , idt(p = -TkVkV, (139) 

and the Lorentz generators are given by 



Y^fei , ^ Mm^ s^Ok 



1 

-^Tk. 



(140) 



In contrast, the multiplications with 7*^ cannot be repre- 
sented by a multiplication of if with a complex 2x2 matrix, 
since 



7 g = (y3 



(141) 



If we express ip* in terms of the two-component complex 
vector 



X = E(p* 



~lT2ip 



-qs + m 
qi - iq2 



the transformation q — > 7*^9 corresponds to 



(142) 



(143) 



We may now introduce the complex four component vector 
*M=(^), (144) 

for which all matrix multiplications q — >■ j'^q can be rep- 
resented by multiplication with complex matrices of the 
Clifford algebra. 



0, 

-i, 



Tk, 

0, 



, 7 = 


-Tk 



^ITk 





0, 

iTk, 
^ki ^ _^_^k 



(145) 



0, T„ 



We can also define the matrix 

7 



-i7°7^7^7^ 



1 



(146) 



for which ip and x ^-re eigenvectors with eigenvalues ±1. 
(Often 7 is denoted as 7^^.) The transformation q ^ Iq 
acts on ipM ^stpM — > il'4'M- Thus the representation (|146p . 
7 = — i/, is consistent with the complex structure. 

Since x is not independent of Lp the spinor ^' m obeys the 
Majorana constraint [loj . 



B = B-'^ = -7^ = 



1 
0-10 
0-100 
10 



(147) 

where Bj^B ^ = (7^)*. The real Dirac matrices 7('jv/) i^- 
the "Majorana representation" © and the complex Dirac 
matrices 7('^^) in the "Weyl representation" (|145l) are re- 
lated by a similarity transformation. 



with 



\/2 



— /l-l-^A" A 



(148) 



/ 1, 0, ,0, 1 



0, 

V 0, 



0, 
1, 




(149) 



and 



A- 




(150) 
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4. Weyl particles 

The propagating degrees of freedom correspond to the 
solutions of the evolution equation. They are most simply 
discussed in terms of the complex equation (|139p . We can 
perform a Fourier transform 

^{t,x) ^ J^0{t,p)e^P^ = J ^0{t,p)e^P^, (151) 

such that the evolution equation becomes diagonal in mo- 
mentum space, 

idtip{t,p) = -pkTk(pit,p). (152) 

The general solution of eq. (|139p obeys 

(f{t,p) = exp{ipkTkt)(f{p), (153) 

with arbitrary complex two-component vectors ifiip). 

We may consider modes with a fixed momentum p. (As 
familiar in quantum mechanics, these modes are not nor- 
malizable for infinite volume and may be considered as lim- 
iting cases of normalizable wave packets.) For every given 
p ^ we can define 

H{p) = ^ , H\p) = 1 , pkPk = \P\', (154) 

such that H{p) has eigenvalues ±1. Decomposing <p{p) 
according to the eigenvalues of H{p), 

H{p)^±{p)=TV±{p) (155) 

we obtain 

<f{t,p) = exp ( - iuj+{p)t)(p+{p) + exp ( - iuj^{p)t)(p^(p), 

(156) 

with dispersion relation 

Lo±{p) = ±\p\. (157) 

We observe positive and negative energies u>±{p). For 
the complex conjugate of ip{x) or (p(j>) the sign of to is re- 
versed. In a four-component representation (jl44p we can 
keep modes with positive energy a; > as the independent 
modes, and relate the other components to them by eq. 
p47p . Thus (fi accounts for one mode ip^ with negative he- 
licity, H(p^ — and x has a propagating mode x_ with 

positive helicity, i?X- = X-- We may associate if with a 
left-handed neutrino or electron, and x with a right-handed 
anti-neutrino or positron. The matrix 7 (|146p projects on 
the left- and right-handed components, 

^(1+7)*A/ , X= ^(1-7)*A/. (158) 

We observe that >p{p) involves for every t two complex or 
four real functions of momentum. This accounts for two 
independent propagating charged particle states. The spe- 
cific state of the propagating modes is determined by the 
initial data >f+{to) and (p_(io) (or equivalently X-(io))- 



We conclude that a particular propagating one-particle 
state can be characterized by the two complex functions 
(fi+ip) and (fi-ip), corresponding to positive and negative 
helicity of the particle. Then ip(t,x) and qy(t,x) can be 
computed according to eqs. (I15ip . (|136p . This fixes the 
Grassmann element gi{t) by eq. (jl26p and therefore the 
coefficients qrit). For the vacuum go = |0) the classical 
wave function qrit) differs from zero only for those bit 
chains r = [n^(a;)] for which precisely one bit is occupied. 
The probabilities of the corresponding classical ensemble, 
Pr{t) = ql{t), are given for each t for which the occupied 
bit is located at position x and of type 7 by the simple 
relation Prit) — q^^t^x). This gives a direct realization of 
a one-particle state in a relativistic quantum field theory 
by a classical statistical ensemble with time evolution ([5^)) . 

As a concrete example we consider <p~{p) nonvanishing 
only for a single momentum p = (0, 0,p3),p3 > 0. This 
implies 

(^i(t,p) ~exp{ip3iM(p- (0,0,p3)) , Mt,p)=0, (159) 
and therefore 

ifiiit, x) = exp{zp3(t + a;3)}, (160) 

and 

qi{t,x) = L-^/^cos{ps{t + X3)}, 
q2{t,x) = L-^/'^sin{p:i{t + Xs.)}, 
q3{t,x) = qi{t,x)^0. (161) 

The probabilities p{x, 7) for the Ising model are given by 

p(x,l) = L-^cos'^{p3{t + X3)}, 

p{x,2) = L-'sin^psit + xs)}, (162) 

and obey 

p{x, 1) + p{x, 2) = L-^ , p{x, 3) = p{x, 4) = 0. (163) 

With equal probability a bit is occupied for every point a;, 
and only the relative probabilities for the occupation of a 
bit of species 7 = 1 or 7 = 2 oscillates in time and in the 
coordinate X3. A state with momentum (pi,0, 0), pi > 0, 
is realized by 

ifiit^x) = Mt^x) = -^L-^/^ exp{ipi{t + xi)} (164) 
v2 

or 

p{x,l) = p{x,3) ^^L-^cos^{pi{t + xi)}, 
p{x,2) ^ p{xA)^\L-^sm^{Piit + xi)}. (165) 

5. Multi-fermion states 

States with arbitrary n describe systems of n fermions. 
This is not surprising in view of our translation of the clas- 
sical statistical ensemble to a Grassmann functional inte- 
gral. For go = |0) there are only n-particle states, and no 
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hole states. They can be constructed by applying n cre- 
ation operators al^{x) on the vacuum. For example, the 
two-fermion state obeys 

= 4s / Qie{t,x,y)a\{x)al{y)gQ. (166) 

V ^ x.y 



Due to the anticommutation relation 

d d 



{a\{x) , al(y)] 



(167) 



the two-particle wave function is antisymmetric, as appro- 
priate for fermions 



qyeit,x,y) = -qe^{t,y,x). 



(168) 



The Grassmann elements g2 describe states with two Weyl 
spinors. A particular class of states can be constructed 
as products of appropriately normalized one particle states 

q^,{x,y) = qW(x)q(2)(y) _ q(l)(y)g(2)(a;). (169) 

General two particle states are superpositions of such 
states. 

It is instructive to specify the probabilities for the Ising 
model that correspond to one particle with momentum 
(0,0, ps) and another with momentum (pi,0, 0). The 
probabilities Pr are now nonvanishing precisely when two 
bits are occupied, and may be denoted by p{x,j,y,S) = 
p{y,S,x,j) = p^s{x,y) for the state where the occupied 
bits are (x, 7) and (y, 5). For the single particle wave func- 
tions ([TFT|) . ([T64| cq. (fT69l) yields 

p(x, 1,2/, 1) = ^L-^ [cos^ {p^it + x^)} cos^{pi{t + yi)} 



+ cos^{pi{t + xi)} cos^{p3{t + ys)} 
-2cos{p3(t + x^)} cos{pi{t + Xi)} 
X cos{p3 {t + y3)}cos{pi{t + yi)}]. 



(170) 



The first two terms are not distinguished and represent 
the product of the probabilities p{x, 1) (I162p and p{y, 1) 
(jl65p written in a symmetrized form. The last term in eq. 
(|170p reflects the characteristic interference for quantum 
particles, with a negative sign as appropriate for fermions. 
Similarly, one has 

Pi2{x,y) = ^L-^[cos^{p3{t + X3)}sm'^{pi{t + yi)} 



+ sin^{p3(i + 2/3)} cos^{pi{t + xi} 
-2 cos{p3{t + X3)} cos{pi{t + Xi)} 
X sin{p3(i + y3)} sm{pi{t -I- yi)}] . 



(171) 



while P22{x,y) and p2i{x,y) obtain from pii{x,y) and 
^12(2;, y) by exchanging cos O sin, and we observe 



^ Pjs{x,y) = ^[l-cos{p3(x3-y3)-pi(a;i-yi)}]. 

(172) 



7,(5=1,2 



One also finds 
Pi3{x,y) +P23{x,y) +pi4{x,y) 



-P24{x,y) = (173) 

and ^33(0;, y) = Paix, y) — P34{x, y) — 0. Due to the inter- 
ference term ^ cos in eq. (|172p the sum of the probabili- 
ties over all species for two given positions (x, y),p{x, y) = 
P-fS ix,y) now oscillates in dependence on a; — y around 
a mean value L^^. Summing over positions x or y the in- 
terference term vanishes, p{x) = '^yP{x,y) — L^^, such 
that the probability p{x) to find a particle of an arbitrary 
species at x is uniform. 

One may wonder why the interference term is necessary. 
If we drop the interference term in eqs. (|170p . (jl71l) . (I172p 
we obtain a perfectly valid probability distribution. How- 
ever, it does not obey the simple time evolution law (|87p 
or ([50]) . Interference effects are a necessary consequence of 
this law. It is not obvious at all if interesting and neverthe- 
less reasonably simple evolution laws for Ising type systems 
exist that avoid interference effects. In any case, interfer- 
ence is a very natural consequence of the simple evolution 
law (|90l) . We note that due to interference the probabilities 
of finding a particle of any species at x and another at y 
are correlated. For our example, these correlations persist 
for arbitrary large separations {x — y). If we replace the 
momentum eigenstates by wave packets with finite local ex- 
tension one will find interference and correlations only for 
regions where the two one-particle wave packets overlap. 



VII. SYMMETRIES 

1. Continuous symmetries 

Besides Lorentz symmetry the action ^ is also invariant 
under global 5*0(2) rotations 

■05. = COS a ■01 — sin a V'2 , "02 = sin a 0i -I- cos a ^02, 
03 = cos a 03 — sin a V'4 ; "04 = sin a 03 -I- cos a 04. 

(174) 

This is easily seen from the infinitesimal transformation 

(50^ = a/-y505 (175) 

and the relations P = -\, = [7,Tfc] = 0. These 
rotations carry over, to the one-particle and one-hole wave 
functions and g^. 

Correspondingly, the complex two component wave func- 
tions ip and X transform as 



(176) 



The S0{2) rotations are now realized as U{\) phase rota- 
tions. If we define for a general complex field 77 the charge 
Q by the transformation 



7] — e ^77 



(177) 



we infer that ip carries charge Q — 1, while x opposite 
charge Q = —1. If describes degrees of freedom of an 



16 



electron, x describes the corresponding ones for a positron. 
We note that charge eigenstates exist only in connection 
with a complex structure. The real wave function can 
be encoded both in Lp and x ^nd may therefore describe 
degrees of freedom with opposite charge. 

2. Discrete symmetries 

The action (O or ([TH]), ([21]), (gH) is further invariant with 
respect to discrete symmetries. Among them a parity type 
reflection maps 

ij{x) ^Vi^{x), (178) 

with 

[Pi^ix])^ = {^\sM-^), (179) 

and we note — —1. We will associate V with a CP- 
transformation since it is realized for Majorana spinors 
which we may take as eigenstates of charge conjugation 
C, Ctjjj — tpj . The action ([3]) is further invariant under the 
discrete transformation — >■ iTp. The continuum limit is 
therefore also invariant under a parity type transformation 
where V is replaced by IV. 

The implementation of time reflection is more subtle. It 
is realized by a map that involves Grassmann variables at 
different times. In a discrete setting for time it acts as 

ip{t, x) — > T ip{t, x), 
fip{t, x)= f ip{-t, x) for t even, 
ftPit, x)^- f x) for t odd, (180) 

where the matrix T obeys 

f ^i-f^ , f^f = 1 , f^Tkf = -Tfe. (181) 

The different transformation of spinors at even or odd time 
points for the discrete formulation is crucial, since other- 
wise the invariance of x'^'yi^)'^-yi^ + ^) could not be 
realized by an orthogonal matrix T. The action ((T5)) . (PT|) . 
([M)) is invariant under the time reflection (I180p . The rel- 
ative minus sign in the transformation of the two spinors 
in the action ([3]) - one at even, the other at odd t - has to 
be remembered for the invariance of the continuum action. 
We can combine the transformations (|179p . (|180l) to a PT 
transformation 

i^{t,x) I 'ip{—t,~x) for t even 
V'(<,a;)->- / ^(-i, -a;) for t odd. (182) 

For Majorana spinors this plays the role of CPT symmetry. 

VIII. COMPLEX STRUCTURE 

We have formulated the classical statistical description 
of a quantum field theory for Majorana spinors in terms of 
a real Grassmann algebra. All quantities in the functional 
integral ^ and the action © or (HH]), (gT]), (HJ) are real. 



The introduction of complex Weyl spinors in eq. (I136P or 
(|142p reveals the presence of a complex structure in this 
real formulation. This will be discussed in more detail in 
the present section. 

1. Complex Grassmann variables 

Complex Grassmann variables may be introduced in 
analogy to eq. (|136p 

Cl = -^(Vl + 1^2) , C2 = ^(^3 + l^i)- (183) 

Together with the complex conjugate Grassmann variables 

Cr = ;7|(^l - *^2) , C2 =^(^3~ *V'4) (184) 

we have for every x and t four independent Grassmann 
variables Ci j Ci j C2 , C2 ; which replace i/ji , -02 , "03 : "04 • 

A general element of a complex Grassmann algebra can 
be expanded as 

3c ^ ^ ^ai . . .ctk .ai . . .ai (-^1 ■ • ■ ■^/ci Xi . . . X/) 

k,l 

CaAxi)-.-CaA^k)QA^l)---Q,{xi), (185) 

with complex coefficients c. (For a real Grassmann algebra 
the coefficients c are restricted by g* — g.) From a general 
complex we can obtain elements of a real Grassmann 
algebra as 

1 i 

9 = 2^9c + g*) , g' = -^{gc 9*c)- (186) 

On every factor C the action of / amounts to a multipli- 
cation with i. More precisely, we can interpret eq. (I183P 
as a map ip — ?> CbP] with the property C[^V'] — *C[V']- For 
the complex conjugate one has C*[-fV'] = ~*C*[V']- For the 
infinitesimal transformation (I175P one concludes 

<5[Ca,(xi)...a,(xO] =*aOCc.,(xi)...a,(*0, (187) 

where the charge Q counts the number of factors C minus 
the number of factors ^* for a given term in the expansion 
(|185p . In other words, products with (5+ factors C and Q- 
factors C* are charge eigenstates with Q = Q+ — Q-. States 
with a given Q are degenerate since many different choices 
of Q- lead to the same Q. For a real Grassmann al- 
gebra we use eq. (jl86p . Every term in gc with Q 7^ is 
accompanied by a term with opposite charge —Q in g*. 
Thus the expansion (|185p of eq. (|186l) involves "charge 
eigenspaces" with pairs of opposite axial charge, similar to 
in eq. (|174l) . 

2. Complex Grassmann algebra 

It will often be convenient to use a complex Grassmann 
algebra, where the coefficients c in eq. (|185p are arbitrary 
and charge eigenstates belong to the Grassmann algebra 
also for Q ^ 0. This can be mapped at the end to a real 
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Grassmann algebra by eq. (|186p . The action of Q± on (jc 
is represented as 



with 



and 



d 



J_ / d 



d 



dUy) V2 \dMy) dMy) 



(188) 



(189) 



dCAy) 



dCiiy) 

J_ / d 



d 



V2 \dMy) My) 



(190) 



and similar for d/dC,2{y) and d/dQ{y), with -01.2 replaced 
by -03,4- 

The Grassmann derivatives d/dQa obey the standard an- 
ticomniutation relations 



(191) 



For Q^Q+-Q^ eq. (fT88l) yields 



IjSips{y) 



(192) 



and we recognize the conserved quantity B2 in eq. (I116P 
with 62 = 1. On the other hand, the sum yields the con- 
served quantity Af via Q+ + Q- = B — J\f. 

In the complex basis we can write the time evolution 
equation in the form 



with 



idtgc = Hgc (193) 



^{TkdkCh + ^{r*kdkCh . (194) 



Both (3+ and Q- commute with H. The operators creating 
one hole or one particle states read 



d 



1 

1 



d 



(195) 



3. Similarity transformation 

Defining four component vectors Tp = ("01, "02, "03, "04) 
and Cd = (Ci, C2, Ci*, C2) _= (C, C*) the relations (HH, (HH 
can be written as a similarity transformation 

^ = GO , Cd = = Gt^, 

G^G = 1 , dctG = 1, (196) 



with 



G 



G-i 



1 

1 
71 



/I, 0, 1, 

~i, 0, i, 

0, 1, 0, 1 

V 0, -i, 0, i, 

/I, ^, 0, 

0, 0, 1, i 

1, -i, 0, 
VO, 0, 1, -I, 



(197) 



The Grassmann derivatives d/dCd obey the standard anti- 
commutation relations (cf. eq. (I19ip ) 



d 



dCd,u 



Cd,; 



(198) 



Grassmann bilinears transforms as ip'^Mip — Cj-^cCd with 
Mq ~ G'^MG. In particular, one finds 

G^G = 

such that the action ([3]) reads 



0,1 
1,0 



, {T,)^^G^T,G^ ( ^Jll ], (199) 



C^(5t-rfc9fc)C. 



(200) 



Up to the factor 2, which may be removed by a rescaling 
of C, this is the action for a free Weyl spinor. 

Transferring a linear transformation ^' — to Cd,C'd — 
L'(d, requires the transformation property L' — G^^LG = 
G'^LG. We note the difference between the transformation 
property of the kernel of a bilinear M which involves G"^, 
and the transformation of L which involves G~^ = G^. In 
particular, one finds 



GtTfcG 



Tfe, 

0, r* 



(201) 



This difference is connected to the additional factor G^G 
appearing in the transformation of the bilinear 

{G^G){G^TkG) = {Tk)c. (202) 

If we define bilinears in terms of 



c^ = cJ(G^G) = cj(;;J) = (cJ 



(CO, (203) 



i.e. by C^M^O, then has the same transformation prop- 
erty as L. One finds 



G^/G 



i, 
0, -i 



(204) 



and the generators S of the Lorentz group dTSl) are trans- 
formed to 



G^EG 



E', 
0, E'* 



(205) 
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with 2 X 2-matrices S' given by eq. (|140p . 

More exphcitly, the infinitesimal Lorentz transforma- 
tions of the complex two-component spinors S( = 
— ^emn5^™"C: ^rc represented by the complex 2x2 ma- 
trices S™" given by eq. (I140p . We also observe that 
( = E( ~ —iT2C transforms as d( — iem„(S™")-^C such 
that 

^C^C = C1C2 = I [Ms - Mi + i{Mi + Ms)] (206) 

is a Lorentz scalar. The same holds for CiQ such that 
the bilinears ^li^^ — '02V'4 ^^'^ ^i'04 + '02V'3 ^I's separately 
Lorentz scalars. Also the product 

ClC2CrC2 = ^1^2^3V'4 = ^^e^^'''^o.^0^-,^S (207) 

is a Lorentz scalar. The functional measure obeys 
/ dTp4dip3dip2d'4'i =]^rf'07 

= J dC2dCi j dQdCi = J dCdC. (208) 

4. Complex structure for real Grassmann algebra 

In a real even-dimensional vector space a complex struc- 
ture is given by the existence of an involution K, together 
with a map /, obeying 

= 1 , /2 = -1 , {K,I} = 0. (209) 

The matrix K has eigenvalues ±1 and we may denote eigen- 
states with positive eigenvalues by vu and those with neg- 
ative ones by vi,Kvii = vr,Kvi = —vj. The matrix I is 
a map between Vfj and vj, implying that the number of 
independent vj^ and vj are equal. We can choose the vj 
such that Ivr = vj, Ivj = —vr and use these properties for 
defining a map from the real vectors v to complex vectors 
c = vb. + ivi — c{v), with the properties c{Kv) = [c(f )] , 
c{Iv) — ic(v). An example is the map (|136p acting on 
V = {97} with 

^ = ^=("0: ra) '^ = ^~ (210) 

Complex structures can also be defined for K and / acting 
only on subspaces. 

We may define a complex structure on the space of the 
Grassmann variables, v=^j, using eq. (j210p . This results 
in the map (^[ip] given by eq. (|183p . with 

where C pl^^ys the role of c. A transformation ijj — > Aip is 
compatible with this complex structure if A obeys 

[A,i]=0. (212) 

In this case A can be represented by complex matrix mul- 
tiplication acting on 

C[A^] ^iCM- (213) 



The matrices Tk commute with / and are therefore com- 
patible with the complex structure 

ClTk^l^TkCm. (214) 

Also / and E^'' (cf. eq. ([T])) are compatible with the 
complex structure, where the action of S'^'^ on is given 
by eq. (|140p . The matrices A obeying eq. (|212p form 
a group, since the product of two matrices A1A2 again 
commutes with /. This product is represented by complex 
matrix multiplication of ^1 and>l2, C[^i^2'0] = ^l^2C['0]■ 
ln contrast, the matrices 7'^ anticommute with /, cf. eq. 
(fT4|) . They are therefore not compatible with the complex 
structure (|210p . (No set of four mutually anticommuting 
complex 2x2 matrices exists.) One finds 

C[7V] = -T-2C. (215) 

Eq. (|210p defines a complex structure for the Grassmann 
elements involving only one factor of 'ijj. One can extend 
this complex structure to Grassmann elements with an ar- 
bitrary number of factors ip, except for the elements |0) 
and |1). The complex conjugation corresponding to eqs. 
(|183p - (|185p reverses the sign of all basis elements gr which 
have an odd number of factors ■02 or ■04. We will not dis- 
cuss this issue further. It should be mentioned, however, 
that many different complex structures can be defined for 
the Grassmann algebra. For example, one could use for the 
involution K the conjugation C. 



IX. PARTICLE-HOLE CONJUGATION 

The time evolution of the probability distribution 
{Prit)} for the classical statistical ensemble is invariant 
under the exchange of occupied and empty bits. For ev- 
ery classical state r = [n-y(a;)] we can define the anti-state 
f = [nj{x)] with nj{x) = 1 — n.y{x). The particle-hole 
conjugation map r — > f transforms occupied bits to empty 
bits and vice versa. The particle number of a state t, given 
by Nr — n^{x) ~ m, corresponds for the anti-state 

to Nf = (1 — nj{x)) — B — m. The particle-hole 

conjugate of a probability distribution {pr} exchanges the 
role of T and f, i.e. C({pt}) = {PtI; with p% — pf- 

1. Particle- hole conjugate wave functions 

On the level of the Grassmann algebra the particle- 
hole conjugation maps Grassmann elements gm ^ ±gB-m- 
Two probability distributions or associated Grassmann el- 
ements which are related by particle-hole conjugation obey 
the same evolution equation. This is expressed by the in- 
variance of the evolution operator K, (llOOp under the ex- 
change of ■01 (x) and d/dip-ylx), and in the identical evolu- 
tion equations for q and q in eq. (|134p . 

The particle-hole conjugation of Grassmann elements 
g{4') can be realized by the operation C which maps ip-y^x) 
factors into factors 1, and inserts r]^ip.y{x) for the places 
(7,x) where no such factor was present. In sect. IIIII we 
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have defined the conjugate Grassmann element gr (no sum 
over r) 



|0) , Cgr = CTrgr , Cgr = (Trgr- 



(216) 



We may define the operation K, g ^ Kg, which multi- 
plies each factor tp-y{x) by 77^, cf. eq. (|210l) . This can be 
used to write 



C — KCph , Cpi^ — 1, 



(217) 



where Cph replaces ipji^) o 1 without additional sign fac- 
tors. (The operation K could be used for defining a com- 
plex conjugation, = 1-) The operation Cph is another 
realization of the particle-hole conjugation of the probabil- 
ity distribution {pr}, since C and Cph only differ by signs 
of the components of the wave function . 

The particle-hole conjugation Cph obeys the relations 

Cph^,{x)Cphg = g^^g , Cphg^^Cphg = M^)g. 

(218) 

For an arbitrary Grassmann operator A we define the 
particle-hole conjugate operator A'^ as 



A — CphACph . 



(219) 



Since we can write ^ as a sequence of suitable opera- 
tors ij-yix) and d/dipj{x) we conclude from eq. (|218p that 
the operator A'' obtains from A by exchanging ip^f{x) ^ 
d/dipj{x). We observe 



ph 



d 



-|0) 



Cphipfix) = ~ip^{x). 



(220) 



The conjugate of the totally empty possible vacuum 
state go — |0) is given by the totally occupied sequence 
ga = |1), represented by the Grassmann element 1, with 
Afgo — Bgo = moffo- Correspondingly, the particle wave 
functions q-yix) for an expansion around |0) are mapped 
by Cph into the hole wave functions q-y (x) for an expansion 
around |1), replacing the creation operators aj^ by annihi- 
lation operators a^. The discussion of hole states is then 
completely parallel to the particle states. We observe that 
the particle-hole conjugate states can also be described in 
the expansion around the vacuum |0). For example, the 
one hole state a-y(a;)|l) obtains by applying B — 1 creation 
operators on |0). The direct construction by applying the 
particle-hole conjugation to a one particle state is a sim- 
ple way of constructing these states. In the following we 
will use a generalized definition of one-hole states as the 
particle-hole conjugate of the one-particle states. These 
states are always contained in the Grassmann elements g. 
For a particle-hole symmetric vacuum go this generalized 
notion coincides with the definition (jl27|) . If Cphgo — .9o 
does not coincide with go we replace eq. (|127p by 



g-i(t) = / q-/{x)a^{x)go. 



(221) 



For J\fgo = mogo,J^go = {B — mo) go one has now 

AA.g_i = (B- Too -1)3-1- (222) 

Two states which are related by particle-hole conjugation 
have the same Lorentz transformation properties. This fol- 
lows from the observation that for a Lorentz invariant go 
also ^0 is Lorentz invariant, together with the transforma- 
tion properties of the creation and annihilation operators. 
In particular, the one-hole wave function (^^{x) transforms 
in the same way as the one-particle wave function (x) . 

We may represent the one-hole wave function by a com- 
plex two-component spinor 



ip{x) 



qi{x) 
qsix) 



ihix) 

iqi{x) 



and use similar to eq. (|142l) 



-f?3 + ^qi 

qi - iq2 



(223) 



(224) 



In the representation by the complex two-component 
spinor x one finds 



Tk 
^ki 



-2S 



Ofe 



-[Tk,Ti] 



(225) 



In this representation the particle-hole conjugation Cph 
maps a left handed particle with charge Q = 1 into a right 
handed hole with charge Q = —I 

We could combine the one-particle and one-hole states 
and define a complex four component spinor "^ph similar 
to eq. mil). 



^ph = 



(226) 



The Dirac matrices acting on '^ph have the representation 
given by eqs. p45p . (|146l) . The Majorana constraint (I147P 
does no longer hold, however, since (p and x are indepen- 
dent. The spinor ^ph describes now a type of Dirac spinor 
which contains degrees of freedom for both particles and 
holes, corresponding to the eight real functions q-f{x) and 
q-yix). We may denote the particle-hole conjugate Dirac 
spinor by 



hp 



ph ■ 



(227) 



As it should be, it obtains from ^p/i by q^ q-y. 

One may be tempted to identify ip with the left-handed 
positron and x with the right-handed electron. However, 
the spinors ip and tp have both Q ^ 1, while x s^nd x carry 
both Q = —1. The upper two and lower two components of 
'I'p/i carry a different charge Q, such that x cannot be iden- 
tified with the right handed electron. (This should carry 
the same charge as the left-handed one.) The Dirac spinor 
(|226p describes a left handed particle and a right handed 
hole with opposite charge, including the complex conju- 
gate states. This amounts to two independent Majorana- 
spinors, or equivalently to two left-handed Weyl spinors 
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ip and (p with identical charge. If we want to describe the 
usual electron and positron we would need p and x to carry 
the same charge. 

2. Particle-hole symmetric vacua 

An interesting situation arises for a vacuum state 
which is invariant under the particle-hole conjugation, 



An example is the linear superposition 
5o = ^(|0) + |l)) 



(228) 



(229) 



and we will discuss other examples below. For a particle- 
hole symmetric vacuum obeying eq. (|228l) both particle 
and hole states are present in the expansion around 50 1 ap- 
plying either creation or annihilation operators. (One has 
to adapt the normalization of eqs. (|126l) . (jl27l) .) The par- 
ticular vacuum state (j229l) has not a fixed particle number, 
since A/'|0) = 0,A/'|1) — B\\). It corresponds to a classical 
probability distribution with probability 1/2 if either all 
Ut ~ \ or if all tIt = 0. 

A perhaps more interesting candidate for a particle-hole 
symmetric vacuum is a half-filled state with N — B/2. A 
possible candidate is given in the complex representation 
(fmi) . (TO by 



= X{{Ci{^)C2{x)} , i^(g,,o + 5;„). (230) 

V 



This state is Lorentz invariant and obeys 



B 



CphQcfl — g*fi , J^9cfi — -^QcM, 



Qg. 



B 



c,0 



(231) 



and correspondingly Cph.gl^o = 9c,o, A/'g*o = 
■j(7*o, Qgcfl — ~t5c.oi such that go is invariant 
under particle-hole conjugation and normalized, — |0)- 
In order to show that go is a static state we use dis- 
crete time steps such that g{t + e) obtains from g{t) by the 
replacement 



(232) 



If gc.o{t) has the form (|230l) one infers from eq. (I232p 

gcAt + ^) = n {a'Ci{x)C2{x) + ^ [C2{x){n),, 

X 

-Cl(a;)(Tfe)2a] [C,a{x + Ak) - Caix - Afe)] 



+ ^^(Tfc)la(T-/)2/3 [Caix + Ak) - Caix - Afe)] 



[Cp{x + Ai) ^ C^ix - Ai)]}. 



(233) 



The combinatorics of this expression looks at first sight 
complicated. However, eq. (|233l) contains precisely B/2 



factors C which have to be distributed on B/A sites x. 
For every x we can have at most two factors of since 
Ciix) — CK-^) = 0- The only nonvanishing contribution 
must therefore have precisely one factor of ^1 and one fac- 
tor of C2 for every x, and we conclude 



gc,o(i + e) = >igcfl{t). 



(234) 



We next show that the proportionality constant k is real, 
K* = K. Indeed, an imaginary part of this coefficient could 
only arise from terms in eq. (I233P which involve an odd 
number of r2-factors. On the other hand, all terms involv- 
ing an odd number of factors for any k must vanish. 
Such terms have an odd number of ^-factors displaced by 
one lattice unit ±Afe, and this is not compatible with hav- 
ing two C-factors for every x. Finally, we have seen that 
the time evolution preserves the normalization of g{t) and 
conclude k — 1. Thus go in eq. (|230p is indeed a static 
state, JCgo = 0. In consequence of eq. (|206p the state g^ is 
Lorentz-invariant and we can define one particle and one 
hole states by applying annihilation and creation opera- 
tors. We notice, however, that g^ is not a state with zero 
axial charge Q. 

The classical ensemble which corresponds to the half- 
filled vacuum (I230p has equal static probabilities for 4^/2 
Ising states, while the other Ising states have zero prob- 
ability. The states with nonzero probability can be char- 
acterized as follows: (i) On every lattice site two bits are 
occupied and two bits are empty, (ii) The occupied bits on 
every site are in one of the allowed combinations of species 
(1, 3), (1, 4), (2, 3), (2, 4). This follows from the property 



Ci{x)Ci{x)go^Ux)C2{x)go^Q, 



(235) 



which forbids the species combinations (1,2) and (3,4). 
(iii) Only if-even states have novanishing probability. The 
X-even states have an even number of occupied bits of the 
species 2 and 4, while K-oAA states have an odd number 
of occupied bits of this type. This follows from the map K 
which maps gc.o ~^ g*o a^nd reverses the sign of all terms 
in g with an odd number of Grassmann variables ^^2 or ip^. 

3. Particle and hole states for the half-filled vacuum 

The normalized one particle state in an expansion 
around go (j230p is given by the Grassmann element 



giit) = / Q7{t,x)al{x){gcfl + glo) 

J X 

Similarly, the one hole state reads 



(236) 



dc*ixy 



= / 97(^.2;)aT,(a;)(.9c,o +5c,o) 



(237) 



V2 



{Va{x)Cix)gc,0 + ipUxKa{x)glo}, 



and we infer that ip and ip both carry axial charge Q = 
1 — B/2. For all t one finds the relation 

g_,{t)g,{t) = ^ I {p*Jx)p^{x)+pI{x)p^{x))\Q). (238) 
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For (pa{x) = (pa{x) the one hole state g-i is the particle- 
hole conjugate state of gi, such that with the normalization 
of the wave function 



one has 



i?V<7-iWffiW = i- 



(239) 



(240) 



The one-particle or one-hole states (|236p and (|237p differ 
from the states Gi\0) and t/_i|l). By now we have already 
found four Grassmann elements for which the wave func- 
tion q^{x) or q-yix) obeys the Dirac equation (|135p . They 
are all eigenstates of AA with eigenvalues 1,B — 1, B/2 + 1 
and B/2 — 1. All these states are described by the func- 
tional integral (f29|) for a suitable choice of the boundary 
term g^. They correspond to different time sequences of 
classical probability distributions Prit). The question to 
what extent these different states have to be considered 
as different particles or should be identified with the same 
macroscopic particle is open at this stage. It seems also 
conceivable that many more states obey the one particle 
Dirac equation. We postpone this issue to a future investi- 
gation. At this point we only emphasize that this feature is 
present in the standard quantum field theoretic formulation 
of a Grassmann functional integral for massless Majorana 
spinors. It is not a particularity of the realization by a 
classical ensemble for an Ising type model. 

4. Observables for excitations of the half-filled 
vacuum 

A convenient set of classical observables for the charac- 
terization of excitations of the half filled vacuum is given 
by N^ix) which takes the value {N!^{x))^ = 1/2 if the 
bit (a;, 7) is occupied, and (A'^^(x))^ — —1/2 if it is empty. 
This classical observable corresponds to the Grassmann op- 
erator 



M'{x) 



1 



d 



-^^{x) - ip^{x)- 



d 



2 [9V7(a;) ' ' 'dil:.y{x) 

for which the basis elements g-r are eigenstates 
A/-;(x)g, = {N'^{x))^gr, 



(241) 



(242) 



with eigenvalues ±1/2. (The eigenvalue —1/2 is found if 
g-r contains a factor -0^(2;).) We can write 



M'^{x)^N^[x)-]^. 
A local particle number can be defined as 



Ar'(x)=^AA;(:r) = ^A/-,(x)-2. 



(243) 



(244) 



Its spectrum has the values (-2,-1,0,1,2). We further 
may employ the total particle number 



(245) 



d 


d 

-X- 


dCiix) 


dClix) 


d 


d 


dCiix) 




d 


d 


dC2{x) 


+ dQ{x) 


d 


d 


dC2{x) 


dQ{x) 



with spectrum {-B/2, -B/2 +1,.. .B/2 - I, B/2). 
In the complex basis we find the expressions 

Kix) = \[[ 

and we define the combinations 

N\{x) = N[{x) +N!^{x) = Ni{x) +N2{x) - 1 
d d 

d d 

9^1 (x) dip2{x) 
Nlix) = N^{x) +Ni{x) = N3{x) +U4{x) - 1 
d d 

:C2{x) + —--<:;{x)-i 



iCiix) 


+cr(^)) 


-} 


(^1 (t) 

\^i\X) 


'5I \X)J 


J' 


{C2{X) 


+ Q{x)) 




(C2(X) 


-Q{x)) 








(246) 



with 



dC2ixy^' ' dci{x) 

= 7^7^^3(x) + ^^.Mx) - 1, (247) 

M'{x)=Miix)+Mlix). (248) 

The half filled vacuum state go is an eigenstate of A/'j(a;) 
and 7VJ(a;), 

M{ix)go = , Ml{x)go - (249) 

and obeys therefore also 

M'{x)go = , Af'go = 0. (250) 

On the other hand, go is not an eigenstate of the difference 
observables Af^x) - Af2{x) or ^'3(0:) - AfHx). 

The one-particle state (|236p is an eigenstate of A/"', and 
similarly for the one- hole state (|237p . 

M'gi=gi , M'g-i = -g-i. (251) 

This follows from the commutation relations 

[al{x),Af^{y)] = -al{x)6^,S{x,y), 
[a^{x),Af^{y)] ^ a^{x)6^,6{x,y), (252) 

which imply that at^{x)go and a-y{x)go are eigenstates of 
J^^iy) and A/J(j/), with eigenvalues 1,0 or —1,0, respec- 
tively, 

A/'j(y)at (x)go A^~f5{x,y)a\{x)go, 
^i{y)a\{x)ga ^ Ai~^5{x,y)a\{x)go, 
J\f\{y)a^{x)ga = -A^^5{x,y)a^{x)gQ, 
M^{y)a^{x)go = -Ai^6{x,y)a^{x)go, (253) 
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where 



^T7 = 



A 



1 for 7 = 1,2 

for 7 = 3,4 

for 7 = 1,2 

1 for 7 = 3,4 



(254) 



With 



J\f'{y)a^{x)go = -6{x,y)aj{x)go, (255) 

eq. pSip follows directly. 

The Grassmann operators Af'{x) can be used for defining 
a position operator for a single particle as 

X = 1 xN'{x). (256) 

For one-particle states its expectation value is given by 

{X) = / xw{x) , w{x) = {M'{x)), (257) 

J X 

where according to eq. (|106p one has 

w{x) = j V^~giN'{x)gi (258) 
= / Qe{y)q-f{z) / V'ipgi^^{y)Af'{x)gi^^{z), 

J y,z J 

where 

5i,^(z) = V2a\{z)go (259) 
is normalized according to 

T^i^9iAy)9ini^) = ^ejSiy, z). (260) 

Using eq. (|255p and the normalization (I260p yields 

^(^)=E97(^)=E^^(^)^"(^)- (261) 

7 a 

Thus ^(a;) > and w{x) = 1, and we can associate w{x) 
with the probability to find a particle at position x. 

The probabilistic interpretation of w{x) allows us to de- 
fine operators corresponding to general functions of a;, 

(262) 



f{X) = I f{x)N'{x), 
with expectation value 

{f{X)) = / f{x)w{x). 



(263) 



Since f{X) is a linear combination of operators A/"' (a;) and 
therefore of Af^{x) it has a direct correspondence with a 
classical observable which takes in a state r the value 



{.f{X))^ = f f{x)[Y,{N,{x))^-2]. 



(264) 



It expectation value obeys the standard statistical law 

{f{X))=J2pr{f{X))^. (265) 

T 

We define an operator product between operators 

XiXi = / XkXiM'ix). (266) 



It maps the classical observables characterized by {Xk)T 
and {Xi)t to a new observable specified by {XkXi)r ac- 
cording to eqs. (|266|) . (|264p . In particular, we may define 
the dispersion observable 



A2 = ^(l,-(1,))2 



(267) 



Its expectation value has the familiar probabilistic inter- 
pretation 



/ w{x)^{xk - {Xk)f 

J X , 



(268) 



It seems rather obvious that the product (|266p is appropri- 
ate for the description of the distribution of position of a 
particle. 

The classical observable product reads 

{Xk ■ Xi)r = {Xk)T{Xl)r 

xm{N'{x))^{N'{y))^. (269) 



It coincides with the product (|266p only for those classical 
states r for which 

{N\x))^{N\y))^ = {N\x))^5{x,y). (270) 

The classical states r for which the probability is non- 
vanishing for a pure one-particle state gi are precisely 
of this type. Indeed, the Grassmann basis elements gr 
which contribute to gi,-y{z) in eq. (|259l) obey N'{z)gr — 
gr, M'{y ^ z)gr — 0, such that all involved classical states 
have (iV'(x))^ ^ only for one particular position x = z. 
Extending the definition of the position operator (|256l) to 
more general states the product (I266P will no longer coin- 
cide with the classical product (I269P [6] . 

From eqs. (|257p . (I26ip we infer that the expectation 
value of the position observable can be computed according 
to the standard quantum rule from the wave function in 
position space 

ifiX)) = / J2^Ux)f{x)ipAx). (271) 

■JX Q, 

The product (|266l) corresponds to the operator product in 
quantum mechanics and may therefore be called "quantum 
product" @ . The formalism of quantum mechanics can be 
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used for a representation of eq. (|271|) in an arbitrary basis, 
e.g. the momentum basis where X involves a derivative 
with respect to momentum. 

The definition of the position operator (|256l) has been 
adapted to the one-particle state. For a one-hole state 
we have to replace J^'{x) by —M'{x). For this purpose 
we replace N'{x) in eq. (|256p by the product N{x) = 
N'{x)N' which coincides with N'{x) and —N'{x) for the 
one-particle or one-hole states, respectively. This extends 
the definition of the position observable to states which are 
superpositions of one-particle and one-hole states. We now 
have 

w{x)=Y,{<l^,{x)+e,{x)), (272) 

7 

where the correct normalization of q and q corresponds to 
^^w{x) = 1. 

5. Generalized one- particle states 

The notion of a position observable is appropriate for 
generalized one-particle states. (For generalized two- 
particle states a pair of two positions is appropriate.) For 
such generalized one-particle states one may replace in eq. 
(|256p N{x) or M{x) by some operator J{x) with the prop- 
erties 

/ J(x) = 1 , ( J(x)) > 0. (273) 

J X 

Then w{x) — {J{x)) defines a probability and f{X) can 
be defined according to eq. (I262p . The product (|266l) will 
coincide with the classical product only if eq. (|270p holds 
for J{x). 

There exists a wide class of possible notions of general- 
ized one-particle states and the choice of J{x) is not unique. 
As a first example, which is close to a possible experimental 
setting, we may cover some region of space by a grid of de- 
tectors, with centers placed on a cubic lattice with lattice 
points X and lattice distance a. To any given detector we 
associate a cube Ia{x) covering all points in the intervals 
Xk — a/2 < Xk < Xk + a/2. We define a "detector observ- 
able" Jaix) by (Jaix))^ = 1 if {N'{x))^ ^ for any one of 
the points x € hix), and {Ja{x))^ = if (N'{x))^ = for 
all points in the interval Ia{x). In other words, the detec- 
tor "fires" if a microscopic particle is present at any point 
of its "detection region" Ia{x) 

A classical state r describes a generalized one-particle 
state if the detector observables obey 

Y.{Ja{x))^ = l, (274) 

X 

i.e. if precisely one detector of the array fires. We can 
define the position observable by 

{X,)r = Y,x{Ux))^, (275) 

X 

and use the probability that a detector at x fires 

Wa{x) = {Ja{x)) ^J2pr{Ja{x))^, (276) 



where the sum over r is restricted to the generalized one 
particle states. For the associated Grassmann operators 
Jaix) the basis elements gr are eigenstates with eigenvalues 
1 or according to 

Ja{x)gr= {Ja{x))^gr- (277) 

Another example may describe a microscopic particle in 
a typical situation where many local or "microscopic" par- 
ticles or holes are excited. For this purpose we associate 
to a given state r "function observables" f^^r{x) Q. We 
choose the value of at x to be proportional to the number 
of microscopic particles of species 7 in the interval Ia{x), 

A,.(x) = c. (^7(^)).' (278) 

with normalization Ct chosen such that 

X 7 

We observe that f^^r is negative if more holes than particles 
are present in the interval Ia{x). For a suitable subset of 
states T that we associate with generalized one particle 
states we can define 

{ja{x))^^Y.^lr{xl (280) 

7 

and use eqs. (|275p . (I276P for the definition of a position 
observable. Eq. (I270p is no longer obeyed and the "quan- 
tum product" (I266P for observables differs from the classi- 
cal product (|269l) . If the grid of points x is fine enough we 
can associate f^^r{x) with functions of a continuous vari- 
able X. 

The "mesoscopic wave function" ^(a;) can be viewed 
as a set of classical observables, one for every value of 7 
and X. Derivatives of f^^r with respect to x are therefore 
also classical observables. In particular, the components of 
"mesoscopic momentum" 

{Pk)r ^-Y.Y1 fiA^VisdkfsA^) (281) 

X 7,(5 

are classical observables which take fixed values for a given 
state r of the Ising model. We define a mesoscopic density 
matrix by 

Pis{x,y) = {fy{x)fs{y)) 

= E?'-A,r(^)Ar(y)- (282) 

r 

In terms of p the expectation values of position and mo- 
mentum of mesoscopic particles find an expression familiar 
from quantum mechanics 

(Xk) - TviXkp) , (Pfe) = Tr(Ffep), (283) 

with trace Tr taken in internal space (7, S) and position 
space (x, y). The operators in eq. (|283p are given by 

{Xk)-(s{x,y) = 5^sS{x,y)yk, 

{PkMx,y) = -i^sS{x,y)^. (284) 
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A quantum product between the classical observables 
{Xk)r and (-Pfc)T is induced @ by the product of the asso- 
ciated operators given in eq. (|284|) . This product is non- 
commutative. The expectation values of quantum products 
involving arbitrary powers the mesoscopic observables Xk 
and Pk can be computed in terms of the density matrix 
similar to eq. (|283|) . In contrast, classical products of the 
mesoscopic observables cannot be computed in terms of p. 
We may employ a complex formulation using 

= /l,r(S) + j/2,r(S), 
<^2?r(S) = /3,r(S)+i/4,r(x). (285) 

The complex density matrix 

p^p{x,y)^U^Hx){^'^;\x)y) (286) 

obeys the standard conditions for a quantum density ma- 
trix, i.e. positivity, hermiticity and normalization, 

pt = p , trp = 1. (287) 

In the complex representation we find the standard repre- 
sentation of position and momentum operators of quantum 
mechanics (in units with 1) 

OXk 

and therefore the standard commutation relation between 
Xk and Pk- 

The mesoscopic observables X and P and their quantum 
products are well defined classical observables for the Ising 
type model and the associated Grassmann functional inte- 
gral. The question remains if they can be used in practice 
for a realistic situation of a mesoscopic particle. For ex- 
ample, we may want to describe an atom as a mesoscopic 
particle, with microscopic scale at the Planck length or 
even shorter. Indeed, an atom involves many microscopic 
fluctuations, as gluons in its nucleons etc.. The answer to 
our question depends to a large extent on the time evolu- 
tion of the density matrix (|282p . This is dictated by the 
underlying evolution of the probability distribution {prii)} 
or the classical wave function {qT-it)}, 

T 

f^,r{x,t) = qr{t)f-,A^)- (289) 

For an isolated mesoscopic particle the evolution equation 
dtPjs{x,y) = 2^f^^r{x)fs,T{y)qTKrpqp, (290) 

(cf. eq. (|57)) ) should be described by a von Neumann equa- 
tion where the r.h.s. can be expressed in terms of p. This 
will certainly not be the case for general wave functions 
{qr}- Finding out the characteristic properties of states 
for which the "subsystem" described by p is isolated, i.e. 
for which the time evolution of p can be computed without 



using information beyond p, is presumably less a formal 
task. It rather involves physical understanding as illus- 
trated by our example of atoms. In this context we note 
that there are many different ways 6] for the definition of 
function observables f^ .^(x) which obey the normalization 
condition (I279P and therefore define a density matrix (|282p . 

6. Particle-hole identification 

The wave function for single particle or a single hole obey 
the same Dirac equation (|134p . Similarly, multi-particle 
states and multi-hole states obey the same evolution equa- 
tion. For a particle-hole symmetric vacuum it is there- 
fore easy to describe particle-hole symmetric excitations. 
The eigenstates of the Cp/^-operator are given for the one- 
particle and one-hole states by 

9±,7(a;) = ^(97(2;) ±97(2;))' (291) 

and similarly for multi-particle states. The eigenstates of 
Cph with eigenvalue +1 have the property that the time 
evolution of the Grassmann wave function g{t) is smooth 
for neighboring time points t and t + e, if g{t) contains only 
elements with an even number of Grassmann variables ip2 
or -04. For an odd number it changes sign. This follows 
from eq. (|217p which implies Cgr — fCgr in eq. (I55|) . A 
particle-hole symmetric initial state Cphgin = 5m remains 
particle-hole symmetric for all even times, Cphg{t) — g{t). 

Particle-hole symmetric states are linear superpositions 
of particle and hole wave functions which do not allow 
anymore any distinction between particles and holes. A 
"one-particle state" should therefore be seen as a particle- 
hole superposition. It can be simply obtained by identi- 
fying q-y{x) — q-y{x), with an appropriate normalization. 
In terms of the classical statistical ensemble for the Ising 
type model this simply means that the probabilities for the 
particle and hole excitations of the particle hole symmetric 
vacuum state are equal. 

X. CONCLUSIONS AND DISCUSSION 

In conclusion, the real Grassmann functional integral 
based on the action ([3]), or its regularized discretized ver- 
sion (US]), (pij) . (1^^ . defines a fully consistent quantum 
field theory for free massless Majorana or Weyl spinors. No 
complex numbers are needed at this stage. Nevertheless, 
the time evolution of the associated multi-particle wave 
function is unitary. One-particle wave functions obey the 
Dirac equation. Multi-particle wave functions are totally 
antisymmetric, as appropriate for fermions. 

The central result of this paper states that this system 
can equivalently be described by a classical statistical en- 
semble for an Ising type model, with a suitable time evolu- 
tion law for the probability distribution. This constitutes 
an explicit realization of quantum mechanics by a classi- 
cal statistical ensemble, in contrast to widespread belief 
that this is impossible. When formulated in terms of the 
probabilities the time evolution law is non-linear. It 
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takes, however, a very simple linear form once formulated 
in terms of the classical wave function which corresponds 
to a "root" of the probability distribution. 

We explicitly indicate the probability distribution for the 
Ising model which corresponds to a quantum state for a 
propagating fermion. A particularly simple example is a 
probability distribution that vanishes whenever more than 
one bit is occupied or if all bits are empty. The states 
T of the Ising model for which pr differs from zero can 
then be associated with the single occupied bits (a;,7), 
Pr=p{x,'Y). The normalization of the corresponding prob- 
abilities p{x,^) obeys J2xJ2jPi^Tl) = 1- A state with a 
fixed momentum pa > in the 3-direction is given by a 
probability distribution 

P(0,0,P3)(2;, 1) = L~^COS^{p3{t + X3)}, 

Pio,o,ps){x,'2) = L~^sm''{p3{t + X3)}, (292) 

andp(a;, 3) — p{x, 4) = 0. The probability oscillates in time 
and space between species 1 and 2. For states with two oc- 
cupied bits one finds the interference effects characteristic 
for fermions. 

The evolution of the probability distribution (|292l) can 
be expressed in a simple way by a wave function 

(j(x,l) = L-^^^cos{ps{t + X3)}, 
qix,2) = L-^^^sin{p3{t + X3)}, 

Pixn) = qH^n)- (293) 

(There is a freedom of choice for the overall signs of q{x, 1) 
and q{x,2), while continuity and differentiability of q fix the 
relative signs for different x and t.) The complex structure 
is introduced by a complex wave function 

(fix) = q{x,l) + iq{x,2) 

= L-3/2exp{ip3(t + x3)}, (294) 

with 

^^*(a:)(^(x) = 1. (295) 

X 

It obeys the simple evolution law for a relativistic particle 
propagating in the xa-direction 

dtip = dsip. (296) 

Our setting generalizes this evolution equation for arbitrary 
momentum and spin, as well as for multi-particle states, in 
a Lorentz invariant way. It describes free massless Weyl 
or Majorana fermions. The complex structure of quan- 
tum field theory can be associated in a natural way to the 
equivalence of Majorana and Weyl spinors. 

Our formulation shows the particle-wave duality char- 
acteristic for quantum mechanics. The particle aspects 
correspond to the discrete spectrum of observables. For 
example, the observable N^{x) counts for a given state 
of the classical ensemble if the bit (x,7) is occupied or 
empty. It can only take the values N^{x) = 1,0 - either a 
particle of species 7 is present at x or not. According to 



the standard rules of classical statistics each measurement 
of N.y{x) can only yield the values one or zero. In other 
words, in any given state r this observable has a fixed value 
{N^{x)^_^ — 1, 0. On the other side, the probability distri- 
bution is typically a continuous function if we take the 
continuum limit for time and space. The wave aspects are 
related to the probabilities. The characteristic interference 
phenomena arise for an evolution law that is linear in q. 

In our approach the axioms for the relation between 
the formalism of quantum mechanics and the outcome of 
measurements follow directly from the standard axioms of 
probability theory or classical statistics. All measurements 
should be formulated as measurements of classical observ- 
ables in a classical statistical ensemble. We have shown 
that this can indeed be done for all observables that can 
be constructed as functions of Nj{x). The expectation val- 
ues of all such observables only depend on the probabilities 
Pr, while the signs St of the wave function, qr = St^Jp^^ 
are irrelevant. 

In the continuum limit, functions of the position observ- 
ables can easily be constructed from iV^ (x) . For Ising states 
with only one occupied bit one uses the position observable 

X = ^xN^(x). (297) 

This can be generalized to one-particle excitations of a 
more general vacuum, as the half filled vacuum - in this 
case one may add a constant to N^{x) in order to guaran- 
tee J^-y ^li^) — fo'" the vacuum state. The observables 
f{X) replace in eq. (|297p x by f{x). Expectation values of 
f{X) at arbitrary t can be equivalently computed from the 
classical ensemble for the generalized Ising model or from 
the quantum wave function. The relation between the ex- 
pectation values and the results of measurements can be 
directly inferred from the standard classical statistical set- 
ting. In our setting the interference effects characteristic 
for quantum mechanics do not pose any conceptual chal- 
lenge. For example, the interference pattern of a double slit 
experiment simply reflects the probabilities to find a parti- 
cle at a certain position of the detector. In principle, these 
probabilities can be computed using only the evolution law 
for the probability distribution. In practice, the concept of 
a "classical wave function" {^^(t)}, with {prit)} — {q^{t)}, 
is very convenient, however. 

For sequences of measurements one needs to define the 
appropriate "measurement correlation" which is based on 
conditional probabilities [l|. There are may different con- 
sistent possible choices of measurement correlations or the 
associated observable products. In principle, the choice of 
the correct correlation function depends on the specific set- 
ting for the measurements, reflecting how a first measure- 
ment affects the outcome of a second one. For many situa- 
tions the measurement correlation differs from the classical 
correlation function. We have explicitly indicated an exam- 
ple of a mesoscopic particle where momentum and position 
are both classical statistical observables of the Ising model 
which can be expressed in terms of (x) . For this example 
we have defined a non-commuting product of observables 
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such that the associated correlation functions reflect pre- 
cisely the quantum correlations. 

We emphasize that all formal and mathematical aspects 
of our discussion of the classical statistical ensemble for 
the generalized Ising model are specified by the regular- 
ized quantum field theory for Majorana spinors. Inversely, 
none of our formal findings depends on the map to the clas- 
sical statistical ensemble. They are properties of a standard 
Grassmann functional integral describing a quantum field 
theory for fermions. This concerns, in particular, the is- 
sue of different vacua raised in sect. IIXI Our setting can 
be generalized to Dirac particles in an external electromag- 
netic field which will be described in a separate publication. 
In this case the single particle wave function will evolve in 
presence of a potential according to the Dirac-equation. It 
remains to be seen if the mapping to an equivalent Ising 
type classical statistical ensemble can also be used for prac- 
tical purposes beyond the conceptual advances in the un- 
derstanding of quantum physics. 
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El 2'Kikm 1 ^..^ 
•^XP i oixr I 1 r ^ + IJ'^fc.O mod 2N+1, 



rn=-N 
N-1 



i\ — 1 f 



2N+1 

TTikm 
N 



2NSkfi mod 2N, 



(A.5) 



imply the relations 

^ exp {iAm{q- p)) = L^Sg^p, 

m 



^ cos(Amq) cos(Amp) = —{Sq+pfi +Sq^p), 

m 

sin(Amg) sin(Amp) = -^(-(^g+p.o + Sq^p), 

m 

cos{Amq) sin(Amp) — 0, (A. 6) 



APPENDIX A: DISCRETIZATION 



1. Lattice Fourier transform 

The lattice Fourier transform expands (we omit the in- 
dex 7) 

iprn = i— ^ { cos{ Arnqjipniq) - sin(Amg)?/;/(9)} 



with Fourier modes ipRj^q) constrained by 



(A.l) 



(A.2) 



The number of different jk is L, such that the number of 
Grassmann variables ipRil)-! ''PiIq) ^ot each species remains 
L^, since the doubling into and ipj is compensated by 
the identifications for ip{q) and ip{—q) in eq. (jA.2p . The 
transformation (|A.ip can therefore be regarded as a simi- 
larity transformation among the Grassmann variables. 
We note that qk = oi qk = tt/A is not contained among 
the possible values in eq. (lA.ll) . The minimal and maximal 
values oi Qk are —tt/A + tt/I and n/ A — n/l. In the contin- 
uum limit this amounts to \qk\ < tt/A. We may consider 
Qk as periodic variables by identifying qk and + 27r/A. 
The mode with j ~ (0, 0, 0) corresponds to 
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cos 



Tt{xi + X2+ X3) 



2^3 I V I 

Tr{xi + X2 + X3) 



— sm 



I 



(A.3) 



and is antiperiodic in all coordinates Xk- 

We may employ complex variables tp{q) — -^{V'-r('Z) + 
ii/j] {q) } , with a standard complex discrete Fourier trans- 
form 



\ ^ iArnq 



L ^ 2^e 



'i^iq) , i^i-q) = r{q)- (A.4) 



with Sq^p = 1 for jq k = jp.k mod L, and zero otherwise. 
The identities (jA.Sp . (jA.6[) allow us to establish 



^ ip{x)llj{x) = ^ (fimi^m 
X m 



(A.7) 



q Q> 

In the last line we employ (/> which may be defined by the 
condition q^ > 0. The linear transformation between the 
independent spinors ipR.-y{q>)j^i,j{q>) and ipjix) is there- 
fore a rotation with unit Jacobian, such that the functional 
measure in terms of '>pR.i{q>) and '>jj{x) is the same. (This 
may be seen by taking ^ — ^(t + e) and ip = ipit)-) 

2. Next neighbor interactions 

For neighboring lattice sites, with Ak the unit lattice 
vector in the k direction, |Afe| = A, one obtains the relation 

J2^ix)^ix + Ak) 

X 

^\Y{ cos(Agfc) [ipR{q)il;R{q) + (p/ (<?)■(/'/ (9)] 
r(A(jfc) [ipR{q)i>i{q) - ^i{q)il^R{q)\ | 



sm 



(A.8 



This can be used in order to show that the derivative (|26p 
obeys 



^ip{x)dk'4){x) 



(A.9) 
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We may also define 



1 



5kHx) = + ^k) + ^{x - Afe) - 2^(x)}, (A.IO) 

such that 

Y,^{x)~5k^{x) (A.ll) 

X 

= 2k ^ [cos(Aqfe) - 1] [ipR{q)ijFi{q) + ^ i {q)ij} i {q)] . 
q 

Our version of a lattice derivative will be introduced by 
considering the combination of Grassmann variables -Fly (x) 
similar to eq. ((24)) . 



This can be expressed as a linear transformation acting on 

'il){x) 

F^[^{t + t)-x\ =Y,P',s{x,y)Mt + ^,y), (A.13) 



with 

Pis{x,y) = ^^\^{Tk)js[S{x,y ~ Ak) 

k 

- 6{x,y + Ak)] (A.14) 
+aS.ys [Six, y-Ak) + S{x, y + Ak) - 2S{x, y)] |. 

We want P to be an antisymmetric matrix 

Ps-,{y,x)^-P^s{x,y). (A.15) 

This is achieved for an antisymmetric matrix 5*, i.e. Ssj = 
—S-ys- We will further require the properties 

52 = -1 , {5, Tfe} = , 5^ = -5. (A.16) 

This is obeyed for 5 = ±7° or S = ±17°. 

3. Lattice dispersion relation 



The combination in eq. (jA.12[) is chosen such that 
"lattice doublers" are avoided. Our setting can therefore 
describe a single species of Majorana or Weyl fermions. 
With 

J2i'iit,x)Fi[i^it + &,x] 

X q k 

I sm{Aqk)['ipB.,-/{t,q)iTk)y5ipi,5{t + e,q) 

-V'/,7(^' Q){Tk)iS-4'R.s{t + e, q)] 
+a[l -cos(A9fe)] [tpR.'y{t,q)SjSipR.sit + e,q) 

+ipi,-y{t, q)S^5^i^s{t + e, q)] }, (A.17) 



we can express the Fourier transform of P as a matrix in the 
space of {tpR{q),'4'i{q)), which is diagonal in momentum 
space, 



E 



ij^{t,x)F^[tp{t+e);x] = ^^■ijj'^ {t, q)P^s{q)4>s{t+e, q). 

(A.18) 



Here we use 



and P{q) — ~P'^{q) obeys 



Piq) = -]ry^{%~'"'^^ ' ^ff w ),(A.20) 
A V -TfcSfc , a(l - Ck)S ' ' ^ ' 



Ff{x) =^^{{Tk).ysdkips{x) + aSjsSk-ipsix)}. (A.12) with 



Sk = sin(A(7fc) , Cfc = cos(Agfc). 



(A.21) 



The properties (jA.16p of the matrix S ensure that P^ is 
proportional to the unit matrix, 

P"('i)--^{T.'l (9) + ( E [1 - ^'^ (9)] ) ' } • (A.22) 

k k 

All eigenvalues of P^ are thus negative or zero. For a 7^ 
zero eigenvalues only occur if for all k 



Sk{q) = , Ck{q) = 1. 



(A.23) 



This is realized only for q = Q. This contrasts to the case 
a = where the zeros of P{q) for both gfc = and qk = tt/A 
lead to the familiar doubling problem for lattice fermions. 
No doubling is present for a 7^ 0, since 



p2 oU -— 

VA' ' / A2 



(A.24) 



and similar for other modes where qk = tt/A. Except for 
(f = 0, where P{q) — 0, one finds that all eigenvalues of 
P{q) are purely imaginary. 

The matrix P{q) is closely related to the inverse fermion 
propagator in momentum space and to the dispersion rela- 
tion which relates frequencies or energies w to g (cf. main 
text). The energy eigenvalues obey 



{u:^ + P\q))i^{q)^Q 
or a dispersion relation 

Lo\q) = -P\q). 



(A.25) 



(A.26) 



All squared frequencies cj^ are strictly positive except for 
q — Q where w((f — 0). For small one finds the relativistic 
dispersion relation \uj\ ~ \q\. 

4. Lattice derivative and continuum limit 



Since F^{x) as defined by eq. (jA.12p has all the desired 
properties we define the lattice derivative dk by 



Fy{x) 



k 



{Tk)j5dk4's{x), 



(A.27) 
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and identify 

dk^dk-^J2 ^khnTiT^Jsh- (A.28) 

For definiteness we may take a = 1/2. We observe that this 
lattice derivative is compatible with the cubic lattice sym- 
metry of the space lattice. The object 8k transforms 
proportional to a lattice Laplacian and S", / are invariant, 
such that Fj{x) has the same transformation property as 

The continuum limit is associated to the behavior of P{q) 
close to its zeros. It can be taken as the limit A — ^ for 
a fixed value of q, and we concentrate on momenta where 
P'^{q) (|A.22|) remains finite in this limit. This happens 
for (f' ^ A^'^, where in eq. (IA.22P we can replace Sk = 
Aqk , 1 — Cfe = A'^ql/2. The diagonal elements of P{q) 
vanish in the continuum limit and we observe the linear 
dispersion relation 

= E ( nqk, ) ' = (^-29) 

As we discuss in sect. I VII the continuum time evolution 
equation for a one particle wave function q{x) is given by 
the Dirac equation 

dtqix) = Y,Tkdkq{x) = / P{x,y)q{v). (A.30) 

k ''y 

In momentum space, with q{p) = qi{p)\ , this reads 

dtqip) = P{p)q{p). (A.31) 

In a complex formulation with qdp) — -^{qnip) + iqi{p)) 
eq. (|X3T|) yields 

dtqdp) = i^PkTkqdp)- (A.32) 

k 

5. Symmetries 

We emphasize that our definition of the lattice derivative 
(jA.28[) or the corresponding definition (|A.12I) for is not 
compatible with the complex structure (j210l) . No matrix S 
exists which anticommutes with all Tk and commutes with 
/. (Otherwise we could find a complex 2x2 matrix which 
anticommutes with all three Pauli matrices Tk , which is not 
possible.) We will take for definiteness 

S^-f\ {S,i}^0, (A.33) 

which makes the discrete formulation compatible with the 
parity transformation (|179p and time reversal (|180p . The 
action of the lattice derivative ()A.28[) on the complex 
Grassmann variable C becomes then 

dkC = dkC~'-^Y. ^MmTtT:^T2hC ■ (A.34) 



For a ^ this mixes C and C* • 

This incompatibility of dk with the complex structure is 
closely related to the observation that the "regulator term" 
^ a violates the symmetry of continuous rotations (|174p . 
(|175|) . Indeed, the transformation (|175p implies 

= [ - ^ + F^ (/»] . (A.35) 

From {S', /} = one infers 

IF = iY,{Tkdk + aS~5k)^ 

k 

= Y,{Tkdk-aS~5k)H^F[ii:]. (A.36) 

k 

In our case of free fermions this symmetry violation can be 
neglected in the continuum limit. 

The violation of the t/(l)-symmetry acting on C by the 
regulator term refiects the well known fact that the U{1) 
symmetry for a single Weyl spinor may remain anomalous 
for an extended theory even in the continuum limit. Our 
functional measure is invariant under the transformation 
(|174p , (|175l) . If one could find a regulated action for a single 
Weyl spinor that is invariant under U{1) transformations 
this would constitute an anomaly free realization of this 
symmetry. In our setting, the f7(l)-symmetry is preserved 
only for a = 0. In this case, however, we encounter fermion 
doubling and the model no longer describes a single Weyl 
spinor. This is consistent with the expectations from the 
Nielson-Ninomiya theorem [13]. 

6. Rotation of Grassmann variables 

The quantity F^{x) in eq. (l24l) coincides with F^{x), as 
defined by eq. (IA.12p . only in the limit e — > 0. For finite e 
we use the relation (jA.lSp in a generalized matrix notation 

F = P4> (A.37) 

in order to define 

F^-\l- exp(-eP)] 1/; = + Ofe^). (A.38) 
e 

(This yields a corresponding modification of the lattice 
derivative dk for finite e/A.) The reason for this modi- 
fication is that we can now write 

B = ilj-eF ^ exp(-eF)V' = Ri^- (A.39) 

bmce PT - -P we infer that R is indeed a rotation, 
R^R = 1. 



APPENDIX B: FUNCTIONAL INTEGRAL WITH 
CONJUGATE GRASSMANN VARIABLES 

In this appendix we reformulate the Grassmann func- 
tional integral ([^^ with action ([T5 | -(pn | on a coarse 
grained time lattice with time points i,„,m G Xjtm+i — 



29 



tm = 2e. In sect. IIVI we have distinguished between even 
and odd time points. The coarse grained lattice only re- 
tains the even time points. The Grassmann variables for 
odd time points will be associated with a different set of 
Grassmann variables at even t by defining for t even 

ip-y{t + e,x) — i/j-y{t,x). (B.l) 

Since all bilinears in eq. (jl9p involve one Grassmann vari- 
able at even t and another one at odd t, the action on the 
coarse grained lattice will involve terms with one factor ip 
and one factor tp. This formulation will be close to the 
description of classical probability distributions by Grass- 
mann functional integrals developed in ref. f?] , with ip the 
conjugate Grassmann variable of ip- Certain properties, 
as symmetry transformations between ^ and ip, are better 
visible in this formulation. 

1. Conjugate spinors 

We start with the "trivial theory" by setting i^-^ = in 
eq. (j20p . The action can now be written in the form 

tf-2e 

S = Yl Lit), 

X 

+ipj{t + e, x)'i(jj{t + 2e, x)} 
= ^V>^(t,a;)[V'^(t-h?,x)-V7(t,a:)], (B.2) 

X 

with e = 2e and J2t — X^m = If "^^^^ yields the dy- 
namical term of the action in the formalism of ref. [7] . For 
e — > the continuum version reads S = ^^pjdt^j. Here 
we have chosen tj even. We next add the term involving 
F^, namely 

AL{t) - -e ^ {ip^{t, x){Tk)-,sdkMt. 

X 

x){Tk)^sdki's{t + 2e, x)} 
+ {A-l)Y{ip-f{t,x)ij-y{t,x) +i:^{t,x)ip^{t + 2e,x)] 

X 

-e'Y'4!^{t,x){Tk)'^sdkips{t + f-jx) 

X 

+ {A-l)Y^^{t,x)ij^{t + l,x) + C{t), (B.3) 

X 

with = 1- 3e2/(2A2). The correction term C{t) van- 
ishes in the continuum limit e — ^ 0, where 

S= j {iJ-,dt4>-,-i^-,{nUdu4>s]- (B.4) 

Jt.x 

The conjugate Grassmann variables ip^ obey the same 
Lorentz-transformations as ip^. and the manifestly Lorcntz 
invariant form of S can be written as 

S^ - I i^rd^i^ , 4, = ^^7°. (B.5) 

Jt,x 



We observe the Lorentz invariance of J^'ipj{x)'ipj{x). To- 
gether with the invariance of the functional measure 
J VipVip this guarantees that also the measure J Dip — 
J DipVip exp { J^ipj{x)ipj{x)^ is invariant. Therefore the 

relations between ipj and ip-y are the same in all systems 
related by Lorentz transformations. 

2. Wave function 

On the coarse grained time lattice (with only even t) our 
formulation has the structure of a Grassmann functional 
integral with conjugate Grassmann variables, which has 
been discussed extensively in ref. Q- The construction 
of the classical wave function and the classical probability 
density can be taken over. We will omit the correction 
term C{t) in eq. (jB.3p . (This corresponds to a different 
regularization of the functional integral.) The results in 
the continuum limit should not depend on this. We also 
rename e — >■ e, such that the regularized action reads 

L{t') = J^[iP{t',x)[yj{t' + e,x)~yj{t',x)] 

-ei^{t',x)Tkdkip{t' + e,x)y (B.6) 

(We omit to write the spinor indices 7.) As before, we 
also may place a; on a cubic lattice with lattice distance A 
and use the lattice derivative ([?5]) . In eq. (jB.6p we have 
extended the summation to include t' — tf, with 

L{tf)^-J2^itf,x)ijitf,x). (B.7) 

X 

In consequence, the partition function is now defined by 
the functional integral 

Z = JvipViPgf[iP{tf)]e-^g^n[HUn)]- (B.8) 

Here the functional integral / VipVip is an integration over 
all Grassmann variables ip and ip (for all tin < t' < tf). The 
particular probability distribution corresponding to this 
functional integral will be specified by the Grassmann ele- 
ment gin [ipjitin, x)] , which can be considered as a bound- 
ary term. The element gf[ipj{tf,x)] is related to ^[i/'(t/)] 

by the additional integration over ip{tf) - not present in eq. 
p9|) - together with the additional term (jB.7p , 

/ 'Di^{tf)9f[^{tf)] eyip{Y^{tf,x)ip{tf,x)) = .g/. 

(B.9) 

We have seen in sect. IIV| eq. ([85|l that gf is conjugate to 
g\_ip[tf)'\. In turn ['0(i/)] is conjugate to 5 [^/'(iy)] in the 
formulation with conjugate spinors [7,]. As a consequence. 



30 



it obtains from gf [ip^itf^x)] by replacing each Grassmann 

variable ip'y{tf,x) — >■ ■^^(t/, x), and by a total reordering of 
all Grassmann variables. In turn, gf will be related to gin 
as a solution of a time evolution equation and is therefore 
not an independent quantity. 
For some given time t we split 



J X 

s< = J2 ^(^')' (B.IO) 

tin<t'<t 

S> = J2 ^(*')+ / ^'it,x)^{t,x), 
t<t'<tf "'^ 

such that depends only on ■(/'(^in < t' < t) and ^p{tin < 
t' < t) while Sy involves il>{t <t' <tf) and 'tp{t < t' < tj). 
For any given time t we can specify the state of the system 
by the Grassmann element 



ait) = g[m] = / T^^Kt' < tWit' < t)e 



(B.ll) 

It belongs to the Grassmann algebra which can be con- 
structed from ipjit, x) for a fixed t. The conjugate element 
obeys 

m = g[m] = / T^Ht' > t)Vi^it' > t)gfe-^>. (B.12) 



It obtains from g{t) by the replacement — >■ tp{t) and 
total reordering. The proof of this property requires the 
hermiticity of Sm 7\. The quantity g{t) plays the role of 
the wave function for a quantum state in the Grassmann 
formulation. 

As in the main text, the Grassmann element g{t) can 
be expanded in a complete basis of the Grassmann algebra 
which is constructed from Grassmann variables ip-yltyx) at 
a fixed time t. 



g{t) ^ qr{t)gr- 



(B.13) 



For Lr' space points we have B = AL^ Grassmann variables 
and therefore 2^ independent basis elements gr- Each r 
will be associated with a state of a classical statistical en- 
semble, which is characterized by an ordered chain of B 
bits taking the values 1 (occupied or spin up) or (empty 
or spin down), r = [Ua] = [n-^ix)] . For every ipa appearing 
in gr we take the bit a empty, while the bit a is occupied 
if tpa does not appear in g-r- As an example, we can asso- 
ciate the Grassmann element g^ — Vi'04V'5'07 • • • with the 
bit chain t = [0, 1, 1, 0, 0, 1, . . . ]. In our case, the states r 
of the classical statistical ensemble can be identified with a 
three dimensional array of four "species" of bits in a com- 
puter, or with the possible states of a four-component Ising 
model in three dimensions. In terms of the "classical wave 
function" {qr{t)} the probability distribution {prit)} of 
the classical statistical ensemble is given by 



The positivity of Pr is obvious and the normalization 

^p,(t) -^g,(t)q,(i) =AAp(t) - 1 (B.15) 

T r 

can be written in the formulation of this appendix as 

Afpit) = / Dmm9it), (B.16) 



Dipit) = / Vi}{t)Vip{t)cxp{ / ^}^{t,x)'ipy{t,x)} 



Here we use the definition of the conjugate Grassmann el- 
ement 



g{t) ^ qr{t)gr- 



(B.17) 



In eq. (jB.17[) the Grassman element gr is conjugate to the 
basis element gr, obeying 



j Dipgrgp = 5rp 



(B.18) 



(For the example gr = ipi'4'4'4'5'4'7 ■ ■ ■ one has gr — 
. . .ipjipz^iipi.) If g{t) as defined by eq. (|B.12p coincides 
with the conjugate of g{t) as defined by eq. (IB.lSp . it is 
easy to see that Afp{t) is independent of t. The definitions 
(jB.lip and (|B.12I) imply J\fp{t) = Z. In turn, it is sufficient 
to choose a properly normalized gin with AfpiUn) = 1 in 
order to guarantee Z = 1. 

3. Time evolution 

We have seen in the main text that for the action ((T5| - 
((20)) the time evolution of the real vector {qr{t)} is given 
by a rotation, which guarantees the normalization (|B.15p 
if gin is properly normalized. It is instructive to show the 
"unitary time evolution" also in the formalism of this ap- 
pendix in terms of conjugate Grassmann variables "0. We 
start with 



9{t + e) 



dil;{t)dijj{t)e^p{~L{t)}g{t) (B.19) 



dmS{m-y^it + e))e''^git), 



Pr{t)=ql{t). 



(B.14) 



with K, given by eq. (j9ip . Here the Grassmann integration 
extends over all variables ip^lt^x) and ijjy{t,x) at a given 
time t. Similarly, the Grassmann (5-function extends over 
all ^pj{t,x). We recall that g{t) and exp{€JC}g{t) depends 
on the variables ip{t), while g{t + e) involves ip{t + e). The 
shift in variables is provided by the (5-function, which corre- 
sponds to the first term in L{t) (|B.6p (not involving Tkdk)- 
The infinitesimal formulation with continuous time avoids 
the ordering problems of the operator exp{e/C} by neglect- 
ing corrections ~ e. This holds provided that the wave 
function qr{t) depends on time in a sufficiently smooth way. 
(For details of the construction cf. ref. TJ.) 
For our case we infer from eq. (jB.6l) 

5(i + e) = n [ / d^iit,x) I dip^{t,x) (B.20) 
X exp{il;^{t,x)[il;^{t,x) - B^]} g{t), 
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with depending linearly on ip^l + e) 

B^{t + e,x) ^ Ai^-,{t + e,x)-e{Tk)-tsdki^s{t + e,x). (B.21) 
The integral over if)^ yields a (5-function, such that 

X g[i:^{t,x)] ^g[B^]. (B.22) 

In other words, g(t + e) obtains from g(t) by replacing each 
Grassmann variable ■ipj{t,x) — ?> Bj{t + e,x). As in eq. ([SS)) 
Bry is a rotation of such that grit + e) obtains from 
grit) by a rotation, and similarly for qr{t + e). In the 
formalism with conjugate spinors the time evolution of the 
expectation value is given by 

dt{A) = {q[A,K]q) = J Dijg[A,IC]g, (B.23) 

with 

dtg = -IC^g, (B.24) 

and KP" obeying 

' Di^t^gf = J DijglCf (B.25) 

for arbitrary Grassmann elements /. 

4. Symmetries 

The formalism with conjugate spinors is convenient for 
the discussion of those symmetries where Grassmann vari- 
ables at even and odd time points transform differently. 
We will concentrate on continuous symmetries of the action 
which involve transformations between ip and ip. Group- 
ing the four-component vectors ipj and ip-y into an eight 
component real vector ijj, 



(B.26) 



the action is invariant under infinitesimal transformations 
Srp = evVip (B.27) 
if the 8 X 8 matrix V obeys (with Tk = diag{Tk, Tk)) 

[V, T,] = , 1/^ = ^CVC , C = (^°' (B.28) 

Independent matrices obeying the condition (|B.28|) are 

V ^{i , Ci , N , iW}, (B.29) 
where / — diag{I, I), IC — CI, and 



N 



-1, 
0, 1 



w 



0, -1 

1, 



(B.30) 



The transformations generated by V (jB.29p commute with 
the Lorentz transformations. For V — I we recover eq. 
(|175p . while the antisymmetric matrix V — CI = —{CI)'^ 
describes an abelian 50(2) rotation between ip and ip- 

In analogy to eq. (IB.26P we can group the one-particle 
and one-hole wave functions q and g, or ip and tp, into a 
common wave function $. For the complex four component 
spinor $ 





. fqi- 












\93 - 


1- mj 



$ = 



the first two symmetries are given by 



(B.31) 



(B.32) 



The other two possibilities for V are symmetric, N'^ = 
N , (IW)^ = IW, and therefore do not belong to the 
generators of eight-dimensional rotations. For example, the 
transformation associated to the generator N describes a 
scaling 



tp ^ e '^ip , Ip — e^'ip. 



(B.33) 



The action of N and / does not mix ip and ip- It can 
therefore be defined for the Grassmann basis elements gr 
and by extension for arbitrary Grassmann elements g. For 

n 

a Grassmann element 5r = 11 ''Pjki^k) the infinitesimal 

k=l 

transformation associated to N reads 

Sgr ^ -ncNgr = -{B - n)eNgT, (B.34) 

such that Ng-r = {n — B)gr. We can identify N with 
the (shifted) particle number operator (IllSp . N — Af — B. 
The scaling symmetry (jB.33|) is associated to the conserved 
particle number. 



ip{x) Vi>{x) 



(B.35) 



The action ^ is invariant under a particle-hole conjuga- 
tion transformation which exchanges ip-y{x) and 'ip^{x). 



Cph : ^-yix) O ip-yix). 



(B.36) 



This symmetry commutes with the Lorentz transforma- 
tions. It reflects the invariance of the evolution equation 
for the probability distribution of the classical statistical 
ensemble under the exchange of empty and occupied bits. 
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